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ABSTRACT 

Signal  Processing  for  Antisubmarine  Warfare  is  a  short 
course  in  electrical  signal  processing  fundamentals  and  their 
applications  in  the  field  of  antisubmarine  warfare.   It  con- 
tains an  introduction  to  Fourier  transforms  and  their  proper- 
ties, sampling  and  quantization,  filters  and  bandwidth 
requirements,  random  signals  and  noise,  and  an  introduction  to 
four  types  of  processing  equipment;  the  DELTIC,  energy  detec- 
tors, correlation  detectors,  and  beamformers.   Course  objectives 
are  given  in  terms  of  specific  questions  which  a  person  com- 
pleting the  course  should  be  able  to  answer.   The  course  text 
and  illustrative  material  is  contained  in  the  appendix  to  the 
thesis.   The  course  is  designed  to  be  presented  in  the  Fleet 
to  the  personnel  involved  with  the  operation  and  employment 
of  detection  equipment  to  provide  them  a  better  understanding 
of  the  operations  accomplished  by  their  equipment  and  to  develop 
in  them  a  better  appreciation  of  the  problems  and  limitations 
associated  with  signal  detection  in  the  antisubmarine  warfare 
environment. 
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INTRODUCTION 

Signal  Processing  for  ASW  was  developed  as  an  introduc- 
tory course  intended  for  the  Naval  Officer  who  has  had  the 
traditional  training  in  ASW  hardware  operations  and  ASW  ad- 
ministrative requirements.   However,  it  is  applicable  to  any- 
one involved  in  ASW  -  officer,  enlisted,  or  civilian  -  who 
does  not  have  an  engineering  background. 

There  are  several  fields  that  are  of  particular  interest 
to  the  ASW  specialist:  acoustics,  signal  processing,  environ- 
mental sciences,  and  probability/statistics.   This  course 
addresses  only  the  signal  processing  aspects  of  ASW.   Short 
courses  covering  Liie  other  three  aieas  are  also  under  develop 
ment  at  the  Naval  Postgraduate  School.   Other  theses  relevant 
to  these  areas  are  KING  and  SWIFT,  1975;  DECKER  and  SOLLEN- 
BERGER,  1975;  and  SMITH  and  VICE,  1975.   For  anyone  who  takes 
all  four  of  these  courses,  it  is  recommended  that  the  signal 
processing  be  taken  last,  because  the  information  contained 
in  the  other  three  will  enhance  the  student's  ability  to 
obtain  an  understanding  of  the  capabilities  and  limitations 
of  signal  processing. 


COURSE  OBJECTIVES 

The  overall  objective  of  this  course  is  to  familiarize 
the  student  with  the  principles  of  signal  processing  used  in 
ASW  sensor  systems.   The  student  who  successfully  completes 
this  course  will  be  knowledgable  in  the  specific  areas  which 
follow  to  the  degree  indicated. 

A.  SIGNAL  PROCESSING  THEORY 

Understand  the  theory  of  signal  processing  sufficiently 
well  to  do  the  following. 

1.  Explain  the  relationship  of  the  Fourier  Transform  to 
the  "time  domain"  and  the  "frequency  domain". 

2.  Given  a  Fourier  Transform  pair,  identify  the  opera- 
tions they  express. 

3.  Given  a  square  wave  and  a  rectangular  wave: 

a.  sketch  the  correlation  of  one  with  the  other; 

b.  sketch  the  convolution  of  one  with  the  other; 

c.  transform  the  given  waves  from  the  time  domain 

to  the  frequency  domain,  multiply  them,  and  apply  the  inverse 
transform  to  return  to  the  time  domain; 

d.  compare  the  results  obtained  in  a,  b,  and  c. 

B.  SIGNAL  PROCESSING  FUNDAMENTALS 

Be  sufficiently  knowledgable  of  signal  processing  funda- 
mentals to  do  the  following. 

1.   Show  the  difference  between  a  digital  and  an  analog 
signal . 


2.  Discuss  how  noise  limits  quantization  in  the  sampling 
and  quantizing  of  a  signal. 

3.  Relate  Sampling  Theorem,  Nyquist  Rate,  and  "aliasing" 

4.  Given  a  plot  of  signal  and  noise  as  power  density 
versus  frequency,  show  how  filtering  increases  the  ratio  of 
signal  power  to  noise  power. 

5.  Given  a  plot  of  filter  response  versus  frequency, 
identify  a  high  pass,  a  low  pass,  and  a  bandpass  filter. 

6.  Relate  filter  bandpass,  integration  time,  and 
frequency  resolution. 

7.  Define  the  statistical  characteristics  of  a  random 
signal  and  relate  them  to  the  signal  voltage  (or  current) 
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8.  List  the  assumptions  necessary  to  apply  statistical 
techniques  to  the  processing  of  random  signals. 

9.  Contrast  the  cross  correlation  of  correlated  signals 
with  the  cross  correlation  of  uncorrelated  signals.   Show  how 
the  correlated  signals  are  enhanced. 

10.  Explain  the  three  effects  of  filtering  gaussian 
white  noise.   Discuss  the  relationship  of  these  effects  when 
processing  signal  and  noise. 


C.   DELTIC  (DElay  Time  Compressor) 

Be  able  to  describe  the  DELTIC,  including: 

1.  how  the  DELTIC  functions; 

2.  the  benefits  of  using  a  DELTIC  in  the  processing  of  a 
signal ; 


3.   given  the  specifications  of  a  particular  DELTIC,  cal- 
culate the  integration  time  required  and  the  frequency  reso- 
lution of  the  output. 

D.  ENERGY  DETECTION 

Describe  energy  detection  by  explaining: 

1.  how  the  system  determines  that  a  target  is  present; 

2.  how  the  threshold  setting  affects  the  probability  of 
dedection,  P(d),  and  the  probability  of  false  alarm,  P(fa); 

3.  the  significance  of  the  direct  relationship  of  a 
P(fa)  with  a  given  P(d)  as  represented  by  a  Receiver  Opera- 
ting Characteristics  (ROC)  curve. 

E.  CORRELATION  DETECTION 
Be  able  to: 

1.  describe  the  method  used  to  discriminate  between  the 
returning  echo  and  noise; 

2.  describe  the  effect  of  doppler  on  correlation; 

3.  explain  how  doppler  is  determined; 

4.  describe  how  pulse  length  effects  doppler  resolution; 

5.  describe  how  pulse  length  effects  range  resolution; 

6.  list  the  advantages  of  using  an  FM  slide; 

7.  state  the  advantages  and  disadvantages  of  using 
Pseudo- Random  Noise  pulses  for  search  and  tracking; 

8.  given  a  two-hydrophone  array  being  used  with  a  corre- 
lation detector,  sketch  and  explain  how  the  signal  direction 
is  obtained,  why  a  direction  ambiguity  exists,  and  how  the 
ambiguity  may  be  resolved; 


9.   explain  the  limitations  on  the  type  of  signals  which 
may  be  processed  passively  by  a  correlation  detector,  and 
why  these  limitations  exist. 

F.   BEAM  FORMING 

1.  Describe  the  method  used  by  the  beam  former  to  obtain 
direction. 

2.  Name  the  type  detector  used  with  a  passive  beam 
former. 

3.  Given  a  hydrophone  array  in  a  beam  former  configura- 
tion, sketch  the  output  of  the  beam  former  for  a  signal  in 
the  beam  and  a  signal  out  of  the  beam.   Show  by  comparing  the 
outputs  the  importance  of  setting  the  threshold  to  determine 
a  target  present  in  the  beam. 

THE  STUDENT 

It  is  expected  that  the  student  will  be  currently  in- 
volved in  ASW.   Prior  courses  in  probability,  college  algebra 
and  calculus,  and  basic  electricity  are  helpful  but  not  re- 
quired, since  the  course  is  self-contained. 

THE  INSTRUCTOR 

There  are  three  basic  qualification  needed  by  the  instru- 
tor  for  this  course:  a  background  in  electrical  engineering 
that  includes  the  fundamentals  of  signal  processing;  a  know- 
ledge of  current  practices  in  ASW;  the  ability  to  teach.   The 
Naval  Postgraduate  School  has  faculty  available  with  these 


qualifications,  and  therefore  should  be  considered  as  a  logi- 
cal source.   MPS  graduates  may  also  provide  a  source  of  people 
with  the  requisite  qualifications. 

COURSE  DEVELOPMENT  AND  PRESENTATION 

This  course  was  designed  to  provide  a  thorough  understand- 
ing of  the  principles  involved  in  the  application  of  spectrum 
analysis  to  ASW  signal  processing,  and  the  operation  of  the 
DELTIC,  energy  detectors,  correlation  detectors,  and  beam 
formers.   These  processing  methods  comprise  the  core  around 
which  the  U.S.NAVY's  ASW  sensor  systems  are  built.   In  pro- 
ceeding toward  the  above  goal,  the  course  follows  a  logical 
approach  of  building  from  the  basics  of  Fourier  Transforms 
and  Fourier  Transform  Properties.   Following  those  two  sec- 
tions is  a  section  on  Sampling  and  Quantization  and  a  section 
on  Filters  and  Bandwidth  Requirements.   These  latter  two 
sections  may  be  interchanged,  but  both  are  required  for  the 
section  on  Random  Signal,  Power  Spectral  Density,  and  Noise. 
The  methods  of  processing  are  presented  last;  in  order  they 
are,  DELTIC,  Energy  Detection,  Correlation  Detection,  and 
Beam  Forming.   Where  appropriate,  specific  hardware  has  been 
mentioned,  but  it  is  emphasized  that  the  course  is  primarily 
devoted  to  principles  and  methods  that  may  be  applied  to  var- 
ious systems  employed  on  either  air,  surface,  or  subsurface 
platforms. 

This  course  is  designed  to  be  taught  in  the  Fleet. 
Although  it  requires  no  outside  references,  each  student 


should  be  provided  with  a  complete  copy  of  the  course  materi- 
al for  his  (her)  use  and  retention.   The  course  objectives 
are  very  specific.   They  are  designed  for  use  by  the  instruc- 
tor and  student  as  a  guide  for  measuring  achievement:  they 
indicate  to  the  student  what  he  should  learn;  they  should 
serve  the  instructor  as  the  basis  for  examination  questions. 
The  only  training  aid  required  by  the  instructor  is  a  chalk 
board,  as  students  may  follow  along  on  the  illustrations  in 
the  course  material.   It  will  be  helpful  if  35  mm  slides  of 
all  the  illustrations  are  available.   Overhead  projectors, 
although  common  at  shore  stations,  are  not  common  on  board 
ship,  whereas  a  35  mm  slide  projector  is  normally  on  board 
for  other  training  requirements.   References  to  books  which 
cover  the  subjects  in  greater  deLaxl,  and  which  have  been  ox 
assistance  to  the  authors  in  the  preparation  of  this  course 
are  included  in  the  bibliography. 

If  the  course  is  presented  independently,  no  more  than 
two  hours  per  day  should  be  devoted  to  in-class  presentation. 
The  student  should  plan  to  spend  a  minimum  of  one  hour  of  re- 
view and  study  for  each  in-class  hour.   The  student  must  have 
time  for  the  information  to  become  impressed  in  his  memory 
and  become  familiar  to  him.   If  all  four  courses  are  present- 
ed concurrently  (feasible  in  a  training  environment  due  to 
the  availavility  of  the  students  time  and  approximately  the 
same  time  required  for  each  course) ,  no  more  than  one  in- 
class  hour  per  day  should  be  spent  on  any  one  course.   Regard 
less  of  the  method  used,  the  instructor  should  be  available 
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at  other  than  schcduled-class  time  to  clarify  points  on 

which  the  student  may  have  questions. 

The  course  is  designed  for  twenty  in-class  hours.   The 

in-class  time  requirement  schedule  is  as  follows. 

SECTION  HOURS 

Introduction  to  Fourier  Transforms-  --------  4 

Fourier  Transform  Properties-  -----------  4 

Sampling  and  Quantization-  -------------  1 

Random  Signal,  Power  Spectral  Density,  and  Noise-  -  4 

DELTIC-  ----------------------  1 

Energy  Detection-  -----------------  2 

Correlation  Detection-  --------------  2 

Beam  Forming-  -------------------  1 

Tl-_        —  „^-_  —  -5—    T  ~  ~ ,  •. .  ~   J         ~    ~  „,-,,.-:, „„-l„.„.«-         4.^         4-]^^,+.         „  r\yr  r\ -*>  r\  A         T   ^ 
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a    two-credit  hour,  one  quarter  course  at  NPS.   The  sections 
which  required  the  most  time  to  develop  were  the  two 
sections  concerning  Fourier  transforms  and  the  section 
on  Random  signal,  power  spectral  density,  and  noise. 
In  these  sections  the  emphasis  is  on  development  of  theory 
and  concept  which  proved  to  be  the  most  difficult  to 
explain  clearly  and  in  a  manner  to  preclude  misunderstandinj 
The  other  areas  were  easier  to  address  as  they  lend  them- 
selves much  more  readily  to  a  descriptive  explanation. 
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CONCLUSIONS 

This  course  effectively  covers  the  major  areas  of  signal 
processing  applicable  to  ASW.   Its  development  has  been  di- 
rected toward  the  understanding  of  theories  and  principles 
used  in  current  ASW  signal  processing  methods.   These  theo- 
ries and  principles  will  remain,  although  in  future  proces- 
sing methods  the  relative  importance  of  their  individual 
applications  may  be  altered.   For  example,  the  current  trend 
in  ASW  signal  processing  is  toward  more  utilization  of  digi- 
tal processing  methods.   This  will  increase  the  importance 
of  thoroughly  understanding  sampling  and  quantization. 

An  area  which  remains  to  be  addressed  is  the  use  of 
digital  computers  in  signal  processing  methods.   This  should 
include  an  introduction  to  boolean  algebra,  elementary  logic 
and  set  theory,  the  use  of  truth  tables  and  Karnaugh  maps, 
binary  arithmetic  and  electronic  logic  circuits.   It  is 
recommended  that  such  a  course  be  developed  as  part  of  the 
continuing  education  program  at  the  Naval  Postgraduate  School 
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APPENDIX 

I.  INTRODUCTION  TO  FOURIER  TRANSFORMS 

An  introduction  to  signal  processing  requires,  for 
understanding  of  many  processing  schemes,  a  working  fam- 
iliarity with  Fourier  transforms.   To  gain  an  understanding 
of  what  Fourier  transforms  are,  what  they  can  do,  and 
how  to  use  them,  let  us  first  examine  a  common  type  of 
signal.   The  electrical  power  used  in  the  United  States  is 
60  Hertz  (cycles  per  second)  alternating  current.   That  is, 
the  current  changes  its  direction  of  flow  in  such  a  manner 
as  to  complete  a  cycle  sixty  times  per  second.   If  we  plot 
the  amplitude  of  the  current  in  a  given  direction  as  a 
function  of  time,  the  result  is  a  sinusoid. 


Since  this  is  a  description  of  the  signal  as  a  function 
of  time,  it  is  said  to  be  the  "time  domain"  description. 
If  we  plot  the  signal  as  a  function  of  frequency,  we  get 
the  "frequency  domain"  description. 
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Since  the  signal  consists  of  a  single  frequency,  there 
is  but  one  frequency  represented  on  the  f-domain  plot. 

Let  us  now  look  at  a  different  signal,  say  an  alternat- 
ing signal  of  130  Hz.  The  time  domain  and  frequency  domain 
plots  are  as  follows. 


•a 
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?! 
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180 


Note  the  way  in  which  the  changes  in  the  signal  are  evidenced, 
both  in  the  time  domain  and  the  frequency  domain. 
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Likewise,  two  other  signals,  300  Hz  and  420  Hz 


to 
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So  far  we  have  looked  at  signals  of  a  single  frequency. 
We  are  going  to  look  at  signals  made  up  of  many  frequency 
components,  but  first  let  us  review  the  concept  of  phase. 
PHASE:   The  phase  of  a  sinusoidal  signal  with  respect  to  a 
reference  signal  is  the  relationship  between  corresponding 
parts  of  their  cycles  in  the  time  domain.   For  example,  take 
a  signal  for  reference. 


Comparing  a  signal  to  this  reference  we  can  determine 
the  phase  of  the  new  signal  with  respect  to  the  reference. 
If  the  new  signal  crosses  the  t  axis  at  the  same  time  and 
in  the  same  direction  as  the  reference,  it  is  said  to  be 
"in  phase"  with  the  reference,  or  to  have  a  phase  angle  of 
0° .   If  the  new  signal  crosses  the  t  axis  on  the  down-swing 
when  the  reference  signal  crosses  on  the  up-swing,  as  shown 
here,  it  is  said  to  be  "180cout  of  phase".   The  reason  for 
using  "180°  will  be  shown  later  in  this  development.   For 
the  moment,  however,  just  accept  the  fact  that  a  signal  is 
said  to  go  through  a  phase  angle  of  360°  each  complete  cycle 
The  phase  of  a  signal  can  either  lead  or  lag  the  reference. 
That  is,  the  signal  may  cross  the  t-axis  on  the  down-swing 
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before  or  after  the  reference  does.   The  first  case  would 
be  leading  and  the  second,  lagging  with  respect  to  the  ref- 
erence.  Because  of  this,  a  signal  is  not  generally  described 
as  having  a  phase  angle  greater  than  180°  since  a  180°  lead- 
ing signal  looks  just  the  same  as  a  180°  lagging  signal.  For 
example,  a  signal  which  is  leading  the  reference  by  270° 
looks  the  same  as  one  lagging  the  reference  by  90° ,  and  is 
usually  described  as  the  latter. 


Signals  with  Multiple  Frequency  Components:  Returning  to 
the  signal  with  more  than  one  frequency  component,  let  us 
look  at  the  way  in  which  this  would  be  represented 
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in  the  time  and  frequency  domains.   Take,  for  example,  the 
four  frequencies  we  looked  at  earlier:  60,  180,  300,  and 
420  Hz.   If  we  had  a  signal  which  was  composed  of  equal  am 
plitude  components  of  each  of  these  frequencies,  the  fre- 
quency domain  plot  would  look  like  this. 
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Now  let  us  look  at  the  time  domain  plot  of  this  signal, 
starting  with  just  the  60  Hz  frequency  present  and  then  as 
we  add  to  the  signal  judiciously  chosen  amplitude  and  ^hsse 
components  of  the  other  frequencies.   First,  take  the  60  Hz. 

portion. 
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Now  add  one-third  the  amplitude  of  the  180  Hz  signal  but 
shift  it  180°  out  of  phase  with  the  60  Hz  signal. 
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Note  the  fact  that  the  180  Hz  component  was  added  180°  out 
of  phase  shows  up  on  the  frequency  domain  plot  as  a  negative 
amplitude.   The  reason  for  this  will  become  apparent  shortly 
For  now,  however,  let  us  continue  with  the  other  frequencies 
To  the  signal  we  now  add  one-fifth  the  amplitude  of  the  300 
Hz  component,  in  phase  with  the  60  Hz  signal. 
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Now  add  one-seventh  the  amplitude  of  the  420  Hz  component, 
but  180°  out  of  phase  from  the  60  Hz  signal. 
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If  we  continue  to  add  carefully  selected  frequency  components 
to  the  signal,  the  time  domain  signal  becomes  more  and  more 
like  a  square  wave.   With  20  components,  the  wave  looks  like 


this  . 
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From  this,  one  can  see  that  by  careful  choice  of  component 
frequencies,  and  their  respective  amplitudes  and  phases, 
that  a  square  wave  can  be  generated  from  purely  sinusoidal 
components.   If  an  infinite  number  of  components  were  used, 
the  resulting  waveform  would  be  a  perfect  square  wave.   This 
concept  is  the  basis  on  which  a  large  part  of  our  signal  pro- 
cessing theory  rests.   There  is  not  time  in  this  course  to 
demonstrate  it,  but  it  is  a  basic  theorem  of  Fourier  analy- 
sis that  any  shape  waveform  can  be  created  from  combinations 
of  pure  sinusoidal  signals. 

Complex  Quantities:   In  order  to  simplify  the  mathematics 
involved  later  in  the  course,  we  must  introduce  the  concept 
of  complex  quantities.   In  solving  the  mathematical  equations 
describing  waves,  one  encounters  solutions  which  contain  terms 
multiplied  by  /T.   As  one  can  readily  see,  there  is  no  real 
number  which,  when  multiplied  by  itself,  will  give  a  number 
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which  is  negative.   Thus,  these  quantities  are  known  as 
"imaginary".   They  are,  in  actuality,  artifacts  of  the  math- 
ematics and  are  not  directly  related  to  tangible  quantities, 
however,  they  prove  to  be  of  use   in  describing  the  behavior 
of  wave  phenomena. 

If  we  plot  all  "real"  quantities  on  one  axis,  all  "imag- 
inary" quantities  on  another,  and  time  on  a  third  axis,  we 
get  a  coordinate  system  like  so: 


RlK 


One  of  the  results  of  the  solutions  described  above  is  that 
a  sinusoidal  signal  in  the  real  plane  (the  plane  formed  by 


the  real  and  time  axes 
can  be  generated  by  a 
vector  rotating  with  a 
uniform  angular  veloc- 
ity in  the  complex 
plane  (formed  by  the 
real  and  imaginary 
axes) .   Looking  at 
this  vector  down  the 
time  axis  shows: 


-  ^  -■-.. 


Rim, 


T.*\t,t_«\Uk.Vt.H 


22 


As  this  vector  rotates  it  is  also  moving  down  the  time 
axis,  so  the  path  swept  out  by  the  tip  of  the  vector  de- 
scribes a  helical  path  around  the  time  axis. 
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At  this  point  it  should  be  noted  that  there  is  another 
measure  of  angular  rotation  known  as  the  "radian".   There 
are  2tt  radians  in  360°  .   This  measure  is  used  because  it 
simplifies  the  mathematics  involved  in  solving  wave  equations 


2k 


Returning  to  the  vector  generator,  we  see  that  the  number 
of  times -it  makes  a  complete  360°  (2tt  radians)  sweep  in  one 
second  of  time  is  the  frequency  of  the  signal.   The  rotation- 
al velocity  of  the  vector,  oj,  is  the  number  of  radians  swept 
out  by  the  vector  per  second.   Thus  co  =  2tt£  since  it  sweeps 
out  2tt  radians  each  360°  rotation  and  makes  f  rotations  per 
second.   Thus  to  is  known  as  the  "radian  frequency"  of  the 
signal . 

If  we  look  at  the  projections  of  the  vector  on  the  real 
and  imaginary  planes  as  it  rotates,  we  find  that  the  projec- 
tion on  the  real  plane  is  a  cosine  wave  and  the  projection 
on  the  imaginary  plane  is  a  sine  wave.   Also,  we  see  that 
the  vector  rotates  360°  during  one  cycle  of  the  cosine  wave 
on  the  real  plane.   This  is  where  the  use  of  360°  of  phase 
angle  arose. 
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The  vector  has  components  in  the  real  and  imaginary 

planes.   A  mathematician,  by  the  name  of  Euler,  has  shown 

the  relationship: 

±i  G  . 

e  J      =  cos  8  ±  j  sin  8( 
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where  j  is  taken  to  mean  "in  the  imaginary  plane". 

By  vector  addition,  the  vector  generating  our  signal  is 
the  vector  sum  of  the  real  component,  cos  oat,  and  the  imag- 
inary component,  j  sin  oat. 

Vector  =  cos  oat  +  j  sin  oat 
We  see  that  this  is  of  the  form  of  Euler's  identity  if  0  =  oat 
Thus,  the  mathematical  description  of  the  vector  is  e~->    ,iJ  . 
This  vector  is  known  as  a  "phasor". 

We  have  seen  that  cos  oat  is  the  projection  of  the  vector 
tip  path  on  the  real  plane.   Solving  Euler's  identity  for 
cos  oat,  we  get: 

Cos  oat  =  Real  Part  of  [e±J  wtJ 

=  Re[e±j  wt]. 
Lihe  convennon  or  using  e  -    tor  cos  oar  is  used  in 
EE  because  it  suggests  several  properties  of  the  cosine 
not  apparent  in  the  written  notation  "cosine  oat".] 

What  does  this  mean?   The  ±  sign  indicates  that  there  are 
two  solutions.   This  corresponds  to  two  phasors  rotating  in 
the  +oa  and  -co  directions,  respectively. 


£»*IA&IHMW 
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These  phasors  unwind  in  opposite  directions  along  the 
t-axis . 


R»»u 


Rcm. 


Looking  at  the  projection  of  these  phasors  on  the  rea]  plane, 
we  see  coscot. 
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Solving  Euler's  identity  for  Cos  cot,  we  obtain 

ej  wt  +  e-j  cot 
Cos  cot  =  t 

We  can  see  this,  as  the  two  contributions  of  the  phasors 
always  add  in  the  same  direction  on  the  real  plane,  there- 
fore, the  +  sign  in  the  formula.   Looking  at  the  projection 
of  these  phasors  on  the  imaginary  plane,   we  see  sin  cot. 


Ri«.u 


Solving  Euler's  identity  for  Sin  cot,  we  obtain 

Sin  cot  : 


i  co t      - i  cat 

eJ     -  e  J 


2j 

Because  the  two  projections  are  always  in  opposite  direc- 
tions on  the  imaginary  plane,  we  get  the  —  sign  in  the 
formula . 

We  have  shown  how  a  single  frequency  plots  as  a  single 
line  in  the  frequency  domain.  This  relates  to  a  line  for 
each  phasor.  But  we  have  seen  that  the  mathematical  solu- 
tions establish  two  phasors  for  each  frequency,  one  rotat- 
ing at  -ho  and  one  rotating  at  -w  .  Thus,  in  actuality,  we 
have  been  showing  only  one  side  of  the  spectrum,  the  posi- 
tive frequency  side.   The  full  spectrum  of  cos  cot  is  then: 
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The  lines  are  both  on  the  same  side  of  the  co  axis  because  of 
the  +  sign  in  the  formula 

j  oat     -i  cat 

Cos  cot  =  ~ 

2         « 

The  full  spectrum  of  sin  oat  is, 


1 

-UJ 

< 

>                      4-UJ 

Note  that  the  lines  are  on  opposite  sides  of  the  w  axis  due 
to  the  —  sign  in  the  formula 


Sin  oat  = 


j  Wt        -j  OJt 

eJ     -  e  J 

2i 


Note  also  that,  the  lines  correspond  to  projections  on  the 

real  plane  for  cos  cat  and  the  imaginary  plane  for  sin  cot. 

We  have  previously  shown  that  arbitrary  waveforms  can 
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be  made  from  combinations  of  sines  and  cosines  of  the  proper 
amplitudes,  phases,  and  frequencies.   How  do  we  determine 
which  frequencies  and  what  amplitudes  and  phases  to  use?   In 
effect,  we  correlate  the  waveform  with  various  frequencies. 
We  compare  the  waveform  to  a  given  frequency  and  see  how  well 
they  match.   One  measure  of  this  is  the  point  by  point  pro- 
duct of  the  two  averaged  over  a  period.   This  measure,  known 
as  C   (n    frequency  coefficient)  is  the  following: 


Cn  =  ^L  \    x(t)  e"jajnt  dt 


where  x(t)  is  the  waveform  of  interest,  e  -'Wn   is  the  cosine 

wpvp    hpino     mmnarprl        and     the     intporql      cpnrp?     tn     cum     the    r\rc\- 

duct  over  the  period,  T  ,  of  the  waveform.   This  determines  how 

o 

closely  the  waveform  matches  a  given  frequency.   In  order  to 
find  the  components  of  a  given  waveform,  we  must  compare  the 
waveform  to  all  frequencies.   The  "spectrum"  of  the  waveform 
is  then  the  summation  of  C  ' s  multiplied  by  the  frequency 
component  (eJ  n  )  over  all  frequencies. 

00 

•X(f)  =Vcn  e^nt 
n=  -°° 

Using  this  technique,  to  analyze  a  square  wave  train,  we  get 
the  spectrum  as  shown  on  the  following  page. 
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We  find  that  all  the  C  ' s  equal  zero  except  for  a  few.   These 
turn  out  to  be  harmonics  of  the  lowest  component.   This  is 
because  the  signal  is  compared  to  sines  and  cosines.   If  the 
signal  is  compared  to  other  forms  of  signals,  the  C  's  will 
not  necessarily  be  harmonics.   This  technique  of  analysis  of 

IftTO  yrt-F/v  T**^  C*        t   r  f\   r>         H  n,rO  1    O'^O'i         ►"%"*»■        "Pf-sii-v^TO-v*         Ottn         f  hn  r         V*  O  C         r\  r\  <"•  r>  rr\  /~%         Is-  *r~\  r\itrn 

as  Fourier  Analysis.   The  fact  that  some  of  the  C  's  are 

J  n 

negative  indicates  that . the  frequency  corresponding  to  that 

C   is  to  be  180°  out  of  phase  with  the  positive  sign  C   fre- 
n  v  v  6    n 

quencies . 

Let  us  look  at  some  square  wave  pulse  trains  and  analyze 
their  spectra. 

i 
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If  we  solve  for  the  C  's  and  plot  we  get  the  curve  as  shown 

n       t-  t. 

on  the  following  page. 
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If  the  period  of  the  pulse  train  is  increased  but  the  length 
of  each  individual  pulse  is  kept  constant,  the  following 
curves  are  obtained. 
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Note  that  the  point  where  the  envelope  of  the  C  's  crosses 

zero  remains  the  same,  but  the  number  of  C  's  increases. 

'  n 

That  is,  there  are  more  frequency  components  present  in  this 
signal  than  there  were  in  the  first.   Let  us  now  increase  the 
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period  even  more,  still  keeping  the  pulse  length  constant 


^CTb. 
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We  see  that  the  zero  crossings  of  the  envelope  remain  con- 
stant, but  the  number  of  C  ' s  increases  again.   If  we  were 

n  b 

to  increase  the  period  to  the  point  where  there  was  only  one 
square  pulse  left,  that  is,  T  -*■   <»,  we  would  find  that  the 

number  of  C  's  would  be  infinite,  but  the  zero  crossings  of 

n  '  & 

the  envelope  would  remain  the  same. 
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From  these  examples  we  can  see  that  the  zero  crossings 

interval  "of  the  envelope  is  related  to  the  inverse  of  the 

pulse  width,  and  the  interval  between  C  '  s  is  related  to 

the  inverse  of  the  period  of  the  signal.   This,  and  other 

properties  of  the  Fourier  relations,  will  be  discussed  in 

greater  detail  later  in  the  course. 

The  form  of  the  envelope  of  the  C  's  for  the  square  wave 

r  n  l 

appears  so  frequently  that  it  has  been  given  a  name  of  its 
own.   The  form  is  (sin  x)/x,  which  is  called  "sine"  (pro- 
nounced "sin-cee"). 

Fourier  Transforms:   We  have  shown  that  as  the  period 
of  a  signal  becomes  infinite,  its  number  of  spectral  compo- 
nents also  becomes  infinite.   In  this  case,  the  spectrum  of 
a  signal  can  be  characterized  by  the  formula  of  the  envelope 
of  the  C  ' s.   Techniques  have  been  developed  for  finding  the 

formula  of  the  envelope  without  solving  for  the  C  's  and  one 

r  t>  n 

method  of  doing  this  is  known  as  the  Fourier  Transform,  de- 
fined thusly: 


00 

X(f)  =  (  x(t)  e"j  wt  dt 


Note  that  this  is  quite  similar  to  the  formula  used  for  find- 
ing the  C  's.   Instead  of  giving  us  the  individual  C  '  s  how- 
&       n  &     &  n 

ever,  this  gives  us  the  formula  for  the  envelope  of  the  C  's 
'       °  r  n 

in  the  f-domain. 

If  we  compute  the  transform  of  a  single  square  wave  using 
this  formula,  we  get  the  following  result  for  the  form  of  the 
s  o  1  u  t  i  o  n  : 
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v/r.    sin  x 

X(f)  =  -*-  sine  x 

K    J  x 

This  is  the  same  result  we  get  if  we  compute  the  C  •  s  for 
a  square  wave  of  T  -*■   °°  and  solve  for  the  formula  of  the 
envelope. 

Since  we  can  find  the  spectrum,  or  f-domain  characteris 
tics  of  a  signal,  given  its  time-domain  formula,  it  is  only- 
reasonable  that  we  should  be  able  to  do  the  reverse  opera- 
tion.  This  turns  out  to  be  the  case,  and  the  operation  is 
accomplished  by  using  the  "inverse  transform". 

X(f)  e^  wtdf 
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By  substitution  in  the  f-domain  formula  of  a  signal  and 
then  integrating,   one  ran  pet  the  time  domain  characteris- 

o  o  ;  tw» 

tics  of  the  signal.  The  two  transform  formulas,  forward  and 
inverse,  are  called  a  "transform  pair".   One  finds  the  f-do 
main  characteristics  given  the  t-domain  formula,  and  the 
other  finds  the  t-domain  characteristics  given  the  f-domain 
formula.   Later  we  will  see  that  these  are  powerful  tools 
used  in  signal  processing. 

For  the  moment,  let  us  return  to  the  cosine  wave  which 
we  looked  at  earlier  and  see  what  happens  if  we  introduce 
modulation.   That  is,  superimpose  it  on  a  carrier  frequency 
cosine  wave  similar  to  AM  radio  transmissions. 

A  modulated  cosine  wave  looks  as  pictured  on  the  follow- 
ing page  in  the  time  domain: 


35 


I-  T  if'  w  i) 


where  the  envelope  varies  as  the  cosine  wave  we  are  looking 
at  and  the  individual  peaks  are  at  the  carrier  frequency. 
What  does  this  signal  look  like  in  the  f-domain?   We  know 
what  the  f-domain  plot  of  our  cosine  wave  looks  like. 


-1  ■■  i 
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The  carrier  is  just  a  cosine  wave  of  a  different  frequency 
so  its  plot  looks  like  this. 


X      f- 
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Note  that  both  f -domain  plots  were  centered  about  f  =  0. 
Now,  if  we  put  the  two  together,  i.e.,  modulate  the  carrier 
with  our  cosine  wave,  we  get  the  following: 


(vO  -ft-tt      i>         fc-f.    H       *~ 


We  find  that  the  resulting  signal  is  of  the  form  Cosine 
(f   ±  £) .   The  f-domain  plot  is  the  cosine  wave  plot  cen- 
tered about  the  carrier  frequency  values  instead  of  about 
f  =  0.   Thus,  modulation  in  the  time  domain  is  equivalent 
to  a  translation  of  the  spectrum  in  the  f-domain  from  cen- 
tered about  the  f=  0  axis  to  centered  about  the  carrier  fre 
quency.   Let  us  look  at  a  modulated  square  wave  pulse  train 


We  know  that  the  f-domain  plot  of  a  square  wave  is  a 
sine   function.   We  have  just  seen  that  modulation  of  a  sig 
nal  translates  the  f-domain  plot  of  the  signal  and  centers 
it  about  the  carrier  frequency.   Thus  the  f-domain  plot  of 
a  modulated  square  pulse  train  is  as  shown  on  the  following 
page. 
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where  the  envelope  takes  the  shape  of  the  square  wave  spec- 
trum but  is  centered  about  the  carrier  frequency  instead  of 
about  f  =  0. 
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II.   FOURIER  TRANSFORM  PROPERTIES 


The  properties  of  Fourier  Transforms  are  stated  most  suc- 
cinctly in  the  accompanying  Fourier  Transform  Theorems.   Know- 
ledge of  the  theorems  is  important  in  the  interpretation  of 
spectra  as  the  theorems  express  the  relationship  between 
time-domain  and  frequency-domain  operations. 

In  the  introduction  to  Fourier  Transforms  it  was  shown  that 
the  frequency  spectrum,  X(f),  of  a  signal,  x(t),  was  the  Four- 
ier Transform  of  the  signal.   That  is,  X(f)  =  j*[x(t)]  and  in- 
versely, x(t)=3"  [X(f)].   These  relations  are  called  "transform 
pairs"  and  are  more  conveniently  denoted  by  x  (t)  *-*-X(f )  .   The 
Fourier  Transform  Theorems  describe  the  properties  of  various 
Fourier  Transform  pairs, 

A.   LINEARITY  (OR  SUPERPOSITION)  THEOREM 


If  a  signal,  x. (t)  ,  is  multiplied  by  a  constant,  a, 
(as  would  happen  if  an  amplifier  with  gain=a  were 
placed  in  a  circuit)  and  a  signal,  x2(t),  is  multi- 
plied a  constant,  b,  then  their  respective  line  spec 
tra  will  also  be  multiplied  by  the  constant  a  or  b. 


bx2  (t ) 


axx(t)  +  bx2(t)  «-*aX1(f)  +  bX£(f) 
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Linear  combinations  in  the  time-domain  become  linear  combina 
tions  in  the  frequency-domain. 

B.   TIME  DELAY  THEOREM 


Given  the  transform  pair  x(t)^-»X(f)  if 
the  signal,  x(t),  is  delayed  by  t   sec- 
onds, as  happens  when  the  signal  passes 
through  a  capacitive  or  inductive  element, 
to  become  a  new  signal,  x(t-t  ),  then  the 
spectrum  is  modified  by  a  frequency  de- 


pendent phase  shift  to  become  X(f)e 


jut, 


x(t-t0)f->X(f)e"jwto 


XtK^ 


m> 


< — > 
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The  translation  of  a  signal  in  time  changes  the  phase  of  the 
spectrum  but  does  not  alter  the  complex  amplitude. 

To  illustrate  the  Linearity  and  Time  Delay  Theorems,  con 
sider  a  signal  which  is  a  linear  combination  of  delayed  sig- 
nals: 

y(t)  =  2x(t-t0)  -  4x(t-3t0) 

Its  spectrum  is  easily  written  as: 

Y(f)  =  2X(f)e_:ia)to  -  4X(f)e";j(1,-5to 
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SCALE  CHANGE  THEOREM 

It  has  been  shown  that  a  translation  of  the  time  origin 
may  be  accomplished  using  the  Time  Delay  Theorem.   The  time 
axis  may  also  be  expanded,  compressed,  or  reversed  by  an  opera- 
tion known  as  "scale  change".   If  a  signal,  x(t),  becomes  a 
new  signal,  x(at),  then: 

(1)   x(at)  is  compressed  if  a  is  positive,  greater  than  one; 


x(t)  = 


t:o 


-—  a  >  1   x(at)  = 


— 

t-* 

d--« 

(2)   x(at)  is  expanded  if  a  is  positive,  less  than  one*, 


x(t)  - 


t-* 


t-< 


tae 


I   +. 


(3)   x(at)  is  reversed  in  time  if  a  is  negative. 


x(t)  = 


t--o 


a  =  -1  x(at) 


t:0 


In  the  Fourier  Transform  pair, 

x(at)  <,  >  I1.  X(-), 

| a |  means  the  absolute  value  of  a,  that  is,  the  numerical  value 
only  is  used,  the  sign  is  dropped. 

The  scale  change  theorem  expresses  the  property  of  recip- 
rocal spreading.   If  the  signal  is  compressed  in  time  by  the 
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factor  a,  its  spectrum  is  expanded  in  frequency  by  1/a.   If 

the  signal  is  reversed  in  time  (a  =  -1)  its  spectrum  is 

f 
X(— j-)  =  X(-f),  which  indicates  a  reversal  in  phase. 
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The  expansion  or  compression  in  the  time  axis  of  a  signal 
occurs  most  often  in  the  playback  of  recorded  signals. 

D.   FREQUENCY  TRANSLATION  (MODULATION)  THEOREM 


If  a  signal,  x(t),  is  multiplied  by  e^ct 
its  spectrum,  X(f),  will  be  translated  in 
frequency  by  ±fc.   Mathematically  stated: 

x(t)eja).ct^X(f-fr), 


Xtt) 


X(t)e 


ja>ct 


4-> 


X(e-ti 


< > 


X(f) 
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It  is  more  common  to  multiply  a  time  signal,  x(t) ,  by  the  real 
part  of  e^  c  ,  i.e.,  by  cos  wct.   From  previous  work  with 
Euler's  theorem,  it  is  recalled  that 

x(t)  cos  uct  =  l/2x(t)(e:J*a)ct  +  e'^c*)- 
Taking  this  result  and  applying  it  to  the  Frequency  Translation 
Theorem  yields  the  Frequency  Modulation  Theorem: 

x(t)  cos  u)ct  <->4/2X(f-fc)  +  l/2X(£+fc) 


<— > 


Or,  if  a  signal  x(t)  is  multiplied  by  cos  a)  t ,  its  spectrum 
X(f)  is  shifted  up  and  down  in  frequency  by  an  amount  f c .  This 
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formed  by  multiplying  a  rectangular  pulse  train  by  cos  wct. 
As  noted  then,  and  confirmed  by  the  modulation  theorem,  time 
domain  multiplication  becomes  translation  in  the  frequency 
domain.   As  an  example,. a  rectangular  pulse  of  amplitude  A, 
period  t,  centered  at  t  =  0  and  duration  t  is  x(t)  =  An(t/x) . 


X(t) 


<-A 


The  spectrum  is  recognized  as  a  SINC  function,  or 

X(f)  =  A  SINC  fx. 


U-3 


Now,  multiply  that  single  pulse  by  cos  w-t 


-A--> 


x(t)  =  An(t/i)  cos  to  t  . 


-A-** 


By  the  modulation  theorem,  the  spectrum  is: 

X(f)  =  1/2  A  SINC  (f-fc)i  +  1/2  A  SINC  (f+f  )t, 


r-fc 


A* 


^^ 


f+fc 


This  is  the  result  that  is  used  extensively  in  single- sideband 
and  double- sideband  HF  communications. 
E.   TIME  DERIVATIVE  THEOREM 


t  e    -t-v, 


dt 


2  „       4.-1-^—   ^4:   ^ 


signal  x(t),  the  result  is  multipli 
cation  in  the  spectrum  by  jar. 

^  x(t)« >  ju>X(f). 
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Therefore,  differentiation  enhances  the  high  frequency  compo 
nents  in  the  spectrum. 

F.   TIME  INTEGRATION  THEOREM 


Taking  the  integral 


f  a  signal  x(t  ) 


1 


is  equivalent  to  multiplying  by  -A-  (di 
viding  by  jco)  in  the  frequency  domain: 


W_  no 


X(t')  atv*i  x(f). 
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time 


Multiplication 

by  * 
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in  frequency 


Integration,  then,  suppresses  the  high  frequency  components 
but  the  low  frequency  components  are  unaffected. 
G.   Convolution 

Take  two  signals  in  the  time  domain  t?,  for  example,  x(t') 
and  y(t  ).   Now  "flip"  the  y(t')  signal  to  the  other  side  of 
the  axis  and  it  becomes  y(-t')  as  illustrated  below. 
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In  order  to  make  a  comparison  between  the  two  signals,  multi 
ply  them  together  at  every  point  in  time  t   from  -«  to  +  °°. 
Mathematically,  this  can  be  expressed  as 


z  = 


00 

f  x(t' 

+J  _  OO 


)  y(-t')dt'  , 


In  the  above  case,  anywhere  in  the  t   domain  where  x(t  ) 
exists,  y(-t')  =  0,  and  where  y(-t')  exists,  x(t')  =  0.  There- 
fore, their  product  is  zero. 

Since  the  signals  are  separated  by  time  t,  if  a  time  delay 
of  t  is  introduced  making  y(-t  )  become  y(t-t  ),  the  signals 
will  overlap.   If  t,  the  time  delay,  is  varied  from  -°°  to  +°°, 
the  y(t-t  )  signal  will  "slip"  across  the  signal  x(t  ).   The 
result  is  expressed  as  a  function  of  time,  t: 


/l+oo 

z(t)  =  J   x(t')  y(t-t')dt'  . 
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The  above  integral  is  known  as  the  "convolution  integral"  and 
it  expresses  convolution  mathematically,  i.e.,  the  area  of 
the  product  of  the  two  signals.   A  more  convenient  notation 
for  convolution  is: 

z(t)  =  x(t)  *  y(t)  t 

A  graphical  interpretation  of  convolution,  as  shown  on  the 
following  page,  may  help  to  see  what  happens  to  the  signals. 
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Two  signals 
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The  shaded  area  on  the 
time  plot  of  both  signals 
represents  the  area  of  their 
"common  product"  for  a  partic- 
ular value  of  time  delay,  t. 
The  shaded  area  on  the  z(t) 
plot  represents  the  summation 
of  the  area  of  their  "common 
product"  as  the  "flipped" 
signal,  y(-t),  is  "slipped" 
across  x (t) . 
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The  convolution,  i.e.,  the  "flip"  and  "slip",  of  two  rec- 
tangular pulses  produces  a  trapezoidal  pulse  with  a  base  of 
length  t-j  +  t j'      Convolution  of  two  identical  rectangular 
pulses  produces  a  triangular  pulse  with  base  =  2t  and  height  = 
AZt. 
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2(0  =  X«J  4"  VC^ 
in  this  instance, 


z(t)  =  x(t)  *  x(t) 


Now  that  convolution  in  the  time  domain  has  been  explained, 
how  does  it  relate  to  the  frequency  domain?  There  are  two  con- 
volution theorems  which  relate  the  two  domains. 
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Convolution  in  the  time  domain  trans 
poses  to  multiplication  in  the  fre- 
quency domain. 


x(t)  •  y(t)*-»X(f)  *  Y(f) 

Multiplication  in  the  time  domain 
transposes  to  convolution  in  the  fre 
quency  domain. 


For  example,  look  at  the  now  familiar  square  pulse  shown 
on  the  following  page. 
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H.   Correlation 

In  statistics,  when  it  is  desired  to  know  how  closely 
one  distribution  resembles  another,  the  two  are  overlayed  and 
the  areas  they  have  in  common  are  multiplied.   This  "common 
product"  represents  the  degree  of  similarity,  or  how  well  one 
distribution  correlates  to  the  other.   In  signal  processing  it 
is  often  necessary  to  compare  one  signal  with  the  same  signal 
which  occurs  later  in  time  or  with  another  signal.   Signals 
vary  in  time  in  contrast  to  distributions  which  are  constant 
in  time.   Therefore,  in   order  to  completely  overlay  one  sig- 
nal on  another,  displace  the  signal  in  time  by  x  units,  then 
vary  x  to  "slip"  the  displaced  signal  across  the  one  with 
which  it  is  desired  to  correlate.   With  each  different  x  used, 
the  "common  product"  will  be  a  different  value.   In  order  to 
obtain  one  number  to  represent  the  value  of  correlation,  sum 
the  "common  product"  over  all  values  of  delay,  t,  and  aver- 
age it  by  dividing  by  time  T  (or  multiplying  by  ™) .   This 
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defines  the  "correlation  function"  R(x)  and  is  described  math- 
ematically as: 

«(,>.»:*  j;.ct,x(f„dt. 

This  is  known  as  a  "time  average"  and  is  denoted  by 

R(t)  =  <x(t)  x(t+x)>, 

where  the  <^  y    (brackets)  indicate   *™  i  jT(x(t)  x(t+x)dt. 
If  the  signal  is  periodic  in  TQ  it  need  only  be  averaged  over 
one  period,  TQ ,  instead  of  all  time.   The  correlation  function 
then  becomes: 
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R(t)  =  j-     /T  x(t)  x(t+x)dt  . 

R(t)  is  more  precisely  referred  to  as  the  autocorrelation  func- 
tion, that  is,  the  correlation  of  one  signal  with  itself.   It 
may  be  desirable,  as  mentioned  earlier,  to  correlate  one  sig- 
nal, x(t),  with  another  signal,  y(t),  displaced  t  units  in 
time  to  become  y(t+x).   This  results  in  the  cross  correlation 
function,  R   (t)  =  <^x(t)  yCt+x^S  .   The  autocorrelation  func- 
tion  is  more  commonly  used.   Where  "correlation  function"  is 
used  it  is  assumed  to  be  the  "autocorrelation  function". 

For  a  graphical  representation,  the  familiar  square  pulse 
is  used.   This,  as  shown  on  the  following  page,  represents  the 
result  of  the  summation  of  the  "common  product"  as  the  delayed 
pulse  is  "slipped"  across  the  other  pulse.   This  initial  part 
of  correlation  is  like  the  previously  described  convolution 
except  that  the  signal  is  not  flipped. 
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The  next  step  is  to  take  the  time  average  of  the  summation  of 
the  common  products.   It  is  seen  that,  as  in  convolution,  cor 
relation  increases  the  width  of  the  function  by  T,  +  T~,  or, 
in  the  case  of  a  square  pulse,  from  T  to  2T. 
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Sum  of 
common  products 


To  time  average,  take  the 
area  =  1/2  bh,  or  l/2(2T)b2, 
and  divide  by  time,  2T, 

l/2(2T)b2    b^ 
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What  does  the  correlation  of  a  sinusoidal  signal  look  like? 
t  us  see  on  the  following  page. 
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When  both  signals  overlay 
completely,  the  value  is 
large  and  positive. 
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A  small  shift  will  give 
some  negative  area  but  the 
result  is  still  positive. 
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A  90°  shift  makes  the  ~ 
positive  and  negative  areas 
equal  and  the  resultant/  val 
ue  zero. 


X 


t-r 


A  shift  of  180°  results  in 
the  entire  area  being  nega 
tive . 


rtf-» 


At  270°  the  positive  and 


negative 


equal  and  the  resultant 
value  is  zero. 


At  360°, a  complete  one  cycle 
shift,  the  signal  is  back  in 
phase  and  the  result  is  a 
maximum  positive  value. 
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The  correlation  function  of  a  sinusoidal  signal  is  anoth 


er 


sinusoid,  and,  if  the  time  average  is  taken,  has  an  ampli 


tude  of 
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Up  to  now  there  have  appeared  many  similarities  between 
convolution  and  correlation.   The  graphical  interpretation 
has  also  shown  this  similarity  but  much  of  that  is  due  to  the 
choice  of  symmetrical  signals  for  illustration.   Now,  look  at 
the  side  by  side  convolution  and  correlation  of  a  saw-tooth 
signal . 


Convolution 
"flip"  then  "slip" 


Correlation 
"slip"  only 
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The  final  results  are  indicated  in  the  next  diagram 


When  "time  Averaged", 
becomes : 


So  it  is  seen  that  if  an  asymetrical  signal  is  used,  the  final 
result  obtained  by  convolution  may  be  quite  different  from  the 
final  result  obtained  by  correlation. 
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III.   SAMPLING  AND  QUANTIZATION 

In  the  real  world,  most  of  the  signals  of  interest  to  an 
operator  of  a  detection  system  are  of  a  continuous,  or  analog, 
nature.   On  the  other  hand,  most  of  the  more  sophisticated 
processors  in  use  today  are  digital  in  operation,  and  thus  can 
not  use  analog  inputs  directly.   Some  means  of  digitizing  the 
analog  signals  must  be  used. 

To  digitize  a  signal,  one  must  describe  the  value  of  the 
signal  at  any  given  time  by  a  digital  number.   This  requires  a 
measurement  and  conversion  of  the  signal  value  to  digital  units 
The  question  of  how  accurately  we  must  measure  the  value  of 
the  signal  in  converting  it  to  digital  form  must  be  considered. 
Take  the  foiluwing  signal,  for  example : 


This  signal  has  a  maximum  value  of  2.5  and  a  minimum  val- 
ue of  0.   Since  we  must  convert  the  continuous  values  of  the 
signal  to  a  set  of  discrete  values,  we  must  then  decide  how. 
many  discrete  values  we  need  in  the  range  from  0  to  3  to  de- 
scribe the  signal  adequately.   If  we  decide  to  use  four  levels, 
that  is,  digitize  to  the  nearest  integer  value,  the  signal 
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will  have  the  following  appearance 

3f 


If  we  double  the  number  of  values  and  digitize  to  the  near 
est  half  integer  value,  the  digitized  signal  looks  like  this. 
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Note  that  as  the  number  of  values  used  in  digitizing  the 
signal  increases,  the  digitized  curve  more  and  more  closely 
approximates  the  original  signal.  The  values  used  in  digi- 
tizing the  signal  are  known  as  quantization  levels.  In  the 
case  of  a  noise-free  signal,  one  can  digitize  the  signal  to 
any  degree  of  accuracy  by  increasing  the  number  of  quantiza- 
tion levels. 
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The  case  of  a  signal  with  noise  is  different,  however. 
In  this  case,  there  is  a  random  fluctuation  of  the  signal 
value  because  of  noise  even  when  the  value  represented  by 
the  signal  itself  is  unchanging.   The  range  of  these  noise 
fluctuations  determines  how  accurately  we  can  measure  the 
value  represented  by  the  signal  even  if  we  have  a  perfect 
measuring  device.   To  illustrate  this,  consider  first  a  sig- 
nal without  noise  and  then  one  with  noise. 
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WITHOUT      NOISE 


with  noise 

We  observe  that  in  the  noisy  signal,  the  value  of  the 
signal  fluctuates  about  some  mean,  or  average,  value  within 
some  limits.   If  the  noise  is  truly  random  in  nature,  the 
mean  value  about  which  it  fluctuates  is  then  the  actual 
value  of  the  signal  being  represented.   At  any  given  instant, 
any  value  of  the  signal  between  the  upper  and  lower  limits 
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of  the  noise  fluctuation  could  stand  for  the  actual  value 
being  represented.   The  accuracy  with  which  this  signal 
can  be  measured  is  then  limited  to  the  range  of  fluctuation 
due  to  noise.   For  example,  suppose  the  signal  of  interest 
is  a  constant  3.0  volts,  but  because  of  noise  the  signal 
fluctuates  between  3.5  and  2.5  volts.   The  most  accuracy 
with  which  we  can  measure  this  signal  is  to  the  nearest 
integer  volt. 

Once  we  determine  how  many  quantization  levels  are 
necessary,  we  mjst  then  determine  how  often  we  need  to  sam- 
ple the  signal  in  order  to  reconstruct  it  accrately. 
Obviously,  if  the  signal  never  changes  value,  but  remains 
constant,  only  one  sample  is  needed  to  reconstruct  the  sig- 
nal for  all  times.   If,  however,  the  signal  is  not  constant, 
we  face  a  real  problem. 

How  often  do  we  have  to  sample  a  given  signal?   Well, 

we  know  that  one  sample  is  not  enough  if  the  signal  is  not 

constant,  and  we  also  know  that  if  the  signal  is  sampled 

continuously  (i.e.,  in  an  analog  fashion)  we  can  reproduce 

it  exactly.   To  sample  a  signal  at  extremely  short  intervals 

is  very  costly  and  difficult  to  accomplish.   What  we  are 

looking  for  is  a  compromise,  a  sampling  rate  high  enough 

that  the  signal  can  be  reproduced  from  the  samples  but  not 

any  higher  than  necessary  in  order  to  keep  cost  and  equipment 

complexity  down. 

We  can  see  intuitively  that  there  must  be  some  relationship 
between  the  rate  of  change  of  the  signal  and  the  sampling  rate 
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necessary  to  be  able  to  reproduce  it.   If  the  signal  changes 
very  slowly,  we  need  sample  it  only  infrequently.   If  it 
changes  rapidly,  however,  we  must  sample  it  frequently  enough 
that  it  does  not  have  time  to  make  several  changes  in  the  inter- 
val between  samples  or  we  lose  these  changes  and  are  thus  unable 
to  reproduce  the  signal.   If  the  sample  intervals  are  too  long, 
the  signal  is  said  to  be  "undersampled" .   This  condition,  in 
addition  to  causing  us  to  lose  some  of  the  information  in  the 
signal,  can  also  cause  other  complications.   Perhaps  you  have 
noticed  that  in  some  of  the  old  cowboy  movies  the  stagecoach 
wheels  appear  to  turn  backwards  at  times.    This  is  a  case  of 
undersampling  which  illustrates  one  of  the  other  effects  called 
"aliasing".   To  explain  this,  let  us  lock  at  a  sine  wave. 
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If  we  sample  this  at  a  frequency  less  than  that  necessary  for 
reproduction  of  the  signal,  we  can  get  a  series  of  samples  which 
give  the  impression  of  a  sine  wave  of  lower  frequency. 
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This  is  a  problem  which  must  be  addressed  when  analyzing  un- 
known signals  for  frequency  components  since  aliasing  can  give 
spurious  results  if  the  analyzer  sampling  rate  is  not  high 
enough.   This  will  be  discussed  further  in  a  minute.   In  the 
meantime,  let  us  look  at  the  sampling  process. 

In  the  time  domain,  we  can  visualize  the  sampling  process 
ideally  as  the  multiplication  of  the  signal  by  a  function  known 
as  the  "ideal  sampling  wave",  which  is  identically  equal  to 
zero  except  at  the  sample  times   ^Ts'  ^Ts'  etc*^* 
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The  duration  of  these  spikes,  known  as  "delta  functions",  is 
infinitesimal,  approaching  zero,  while  their  amplitude  approaches 
infinity,  and  they  are  defined  such  that  their  area  (0  x  °°)  is 
equal  to  1.   When  the  signal  of  interest  is  multiplied  by  the 
ideal  sampling  wave,  we  get  the  product  of  their  areas  at  each 
point  T.   Since  the  area  of  each  delta  function  is  1,  and  its 
width  is  0,  the  product  of  areas  is  equal  to  the  value  of  the 
signal  at  that  time.   Thus  the  products  of  the  two  functions 
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Let  us  look  at  the  frequency  domain  to  see  more  clearly 
what  happens  when  a  signal  is  sampled.   Assume  that  the  signal 
to  be  sampled  is  strictly  "bandlimited" ,  that  is,  it  has  no 
frequency  components  outside  of  some  limit. 
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2W  is  the  bandwidth  of  this  signal.   Although  the  mathematics 
is  rather  involved,  it  can  now  be  shown  that  the  ideal  sam- 
pling wave  spectrum  appears  as  in  the  next  diagram. 
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We  know  from  the  modulation  theorem  that  modulation  of  a  func 
tion  causes  its  spectrum  to  be  shifted  and  centered  about  the 
modulation  frequency  and  otherwise  does  not  change  it.  Thus, 
the  spectrum  of  the  sampled  signal  appears  as  follows. 


From  this  we  can  see  that  for  strictly  bandlimited  signals,  if 
£   >  2W,  there  is  no  overlapping  of  the  spectrum  components 
and  so  the  spectrum  of  the  original  signal  is  exactly  repro- 
duced.  If  the  sampled  signal  is  then  passed  through  a  low-pass 
filter  which  has  the  property  of  allowing  only  those  frequencie. 
below  its  cutoff  to  pass,  we  can  exactly  recover  the  original 
signal . 
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It  follows  from  the  above  that  the  sampling  frequency  must  be 
>  2W,  i.e.-,  greater  than,  or  equal  to,  two  times  the  highest 
frequency  component  of  the  signal  if  we  are  able  to  reconstruct 
the  signal  from  the  samples.   If  it  is  not,  the  spectrum  of 
the  sampled  signal  has  overlaps  which  introduce  errors  into  . 
the  reconstruction  of  the  signal. 
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These  overlaps  are  the  cause  of  "aliasing"  and  other  incon- 
sistencies in  the  reconstructed  signal.   This  result  is  formal- 
ized in  the  Sampling  Theorem,  or  theorem  of  uniform  sampling: 
If  a  signal  contains  no  frequency  com- 
ponents for  |f|  >  W,  it  is  completely 
described  by  instantaneous  sample  values 
uniformly  spaced  in  time  with  period 

T   <  1/2W. 
s  -   ' 

This  rate  of  sampling,  f   =  2W,  is  known  as  the  Nyquist  Rate. 
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It  is  the  absolute  minimum  sampling  rate  from  which  a  signal 
can  be  reconstructed  under  ideal  conditions. 

Until  now  we  have  considered  only  ideal  sampling  and  recon- 
struction.  In  practice,  the  ideal  is  rarely,  if  ever,  realized. 
Practical  sampling  differs  from  ideal  sampling  in  three  obvious 
respects  : 

(1)  Practical  sampling  waves  are  not  composed  of  delta 
functions  but  have  a  finite  width. 

(2)  Practical  filters  are  not  ideal. 

(3)  Real  world  signals  normally  sampled  are  not  strictly 
bandlimited . 

In  practice  these  differences  are  of  small  enough  significance 
that  for  all  normally  encountered  cases  only  the  last  is  of 
real  iip.'o o  r t  anc e  . 

Since  most  real  signals  are  not  strictly  bandlimited,  there 
will  be  some  overlap  in  the  sampled  spectrum.   If  the  spectral 
components  outside  some  nominal  value  W  are  negligible,  however, 
the  signal  can  still  be  adequately  described  for  most  purposes 
by  samples  spaced  T   <  1/2W. 
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IV.  FILTERS  AND  BANDWIDTH  REQUIREMENTS 

The  ocean  is  an  extremely  noisy  environment  in  which 
to  search  for  submarines.   Noise  from  waves,  wind,  breaking 
surf,  shipping,  sonic  mammals,  fish  and  crustaceans, 
seismic  activity,  rain,  etc.,  is  present  to  some  degree 
at  all  times.   How  do  we  go  about  sorting  out  these  dif- 
ferent sources  and  detecting  submarines?   One  way  is  to 
use  filters  to  eliminate  as  much  of  the  noise  as  possible. 

In  general,  there  are  three  classes  of  filters  relative 
to  spectral  response:   low-pass,  high-pass,  and  band- 
pass filters.   Ideal  low-pass  filters  allow  frequencies 
below  their  cutoff  frequency  to  pass,  and  screen  out 
all  above  cutoff.   Their  spectral  response  is  shown  in 
the  following  diagram. 
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High-pass  filters  allow  all  frequencies  above  their 
cutoff  to  pass  and  screen  out  all  below. 
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Band-pass  filters  allow  all  frequencies  between  their 
lower  and  upper  cutoff  frequencies  to  pass  and  screen 
out  all  others. 
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Filters  are  characterized  by  their  upper  and  lower 
cutoff  frequencies,  or  their  center  frequency  and  band- 
pass (or  bandwidth),  B.   The  bandwidth  is  defined  as  the 
passband  width  measured  in  positive  frequency  only, 

that  is,  the  difference  f   -  £n  -  B. 
'  u     1 
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Of  course,  no  real  filters  have  the  sharp  cutoff 
characteristics  of  the  ideal  filter,  so  it  is  customary 
to  define  the  bandwidth  of  a  real  filter  as  the  band- 
width between  the  points  v/here  the  response  has  dropped 
3  dB  from  the  maximum.   These  points  are  known  as  the 
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A  signal  which  is  passed  through  a  band-pass  filter 
appears  as  in  the  next  diagram. 
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As  indicated  earlier,  filters  eliminate  as  much  of 
the  noise  as  possible.   How  does  the  filter  actually 
improve  detection  capability?   To  answer  this  question, 
we  examine  the  way  filters  help  improve  detection  in  a 
large  class  of  detectors,  known  as  energy  detectors. 
We  discuss  the  way  these  operate  in  greater  detail  a 
little  later,  but  the  general  principle  of  operation 
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is  that  the  detector  senses  the  total  power  (or  energy) 
present  in  the  input,  and  compares  this  to  a  "threshold" 
value.   If  the  input  exceeds  threshold,  the  detector 
indicates  the  presence  of  a  signal,  and  if  the  input  does 
not  exceed  threshold,  it  indicates  no  signal  present. 

In  order  to  see  better  how  filtering  helps  to  detect 
signals  in  noise,  consider  the  "power  spectral  density" 
plot  of  an  input.   The  power  spectral  density  is  the  power 
in  a  signal  that  is  carried  by  each  frequency  component. 


The  total  power  in  a  signal  then  corresponds  to  the 
integral  of  the  power  spectral  density  over  all  frequencies 
The  power  spectral  density  thus  gives  an  indication  of 
which  frequencies  carry  most  of  the  power  in  a  signal, 
and  the  relative  amount  of  power  at  a  given  frequency 
compared  to  that  at  others.   The  power  in  a  given  freq- 
uency band  is  then  proportional  to  the  area  under  the 
power  spectral  density  curve  between  the  lower  and  upper 
limits  of  the  band.   Comparison  of  the  power  present 
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in  a  signal  to  the  power  of  the  background  noise  present 
indicates "the  detectability  of  the  signal.   Consider 
a  signal  against  a  background  of  "white  noise"  (white 
noise  is  defined  as  having  equal  amounts  of  all  spectral 
components)  to  see  the  gains  from  filtering. 
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We  can  see  that  if  the  detector  examines  the  entire 
frequency  spectrum  and  senses  the  total  power  present 
to  compare  to  a  threshold  value,  the  presence  or  absence 
of  the  signal  does  not  affect  the  input  significantly. 
On  the  other  hand,  if  we  filter  the  input  and  only  allow 
a  narrow  band  of  frequencies  about  the  frequency  of  interest 
as  the  input  to  the  detector,  then  the  presence  of  a  signal 
will  substantially  change  the  amount  of  power  in  the  input. 
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This,  in  effect,  raises  the  signal-to-noise  ratio, 
and  thus  allows  us  to  detect  weaker  signals  than  we  could 
detect  before  filtering.   In  the  case  of  non-white  noise, 
the  improvement  in  signal-to-noise  (S/N)  ratio  can  be 
even  greater. 


powER 


Kat$£ 


Ambient  noise  in  the  ocean  is  an  example  of  non- 
white  noise,  diiu  CaieiUi  filtering  iun  imp  rove  u-ctcctiori. 
probability  immensely. 


$ 


<-c=> 


Energy  detectors  integrate  the  energy  in  the  input 
over  some  time  period,  some  averaging  the  energy  and 
some  not  averaging.   The  value  obtained  in  either  case 
is  compared  to  the  threshold,  which  is  set  at  the  value 
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corresponding  to  ambient  noise  when  the  signal  of  interest 
is  absent.   Therefore,  one  must  determine  how  long  to 
integrate  the  input  before  comparing  to  threshold.   How 
is  this  determination  accomplished?   We  saw  in  the  section 
on  Fourier  transforms  that  there  was  a  relationship 
between  the  period  of  a  signal  and  the  separation  between 
the  frequency  components  that  make  up  the  signal,  namely, 
the  separation  between  components  is  inversely  related 
to  the  period  of  the  signal. 
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Thus  if  we  assume  that  the  integration  time  we  choose 
is  the  period  of  the  input,  and  that  the  input  exactly 
repeats  itself  every  period  ad  infinitum,  then  the  minimum 
resolvable  frequency  difference  between  components  of  the 
input  is  inversely  related  to  the  integration  time.   The 
longer  the  input  is  integrated,  the  more  accurately  the 
frequency  components  present  are  determined.   If  we 
desire  to  resolve  frequencies  1  Hz  apart,  we  must  integrate 
for  at  least  1  second,  for  frequencies  0.1  Hz  apart,  at 
least  10  seconds. 
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Another  consideration  in  detection  system  design 
is  the  fact  that  most  common  filters  in  military  appli- 
cations today  are  composed  of  reactive  elements,  such 
as  RC/LC/RLC  circuits,  crystals,  or  transducer  elements, 
which  have  finite  time  constants  associated  with  them. 
What  does  this  mean?   It  means  that  these  systems  are 
designed  to  resonate  at  a  particular  frequency  or  over 
a  band  of  frequencies,  and  to  block  or  short  out  others. 
These  reactive  elements  require  a  finite  time  to  "ring 
up",  or  establish  resonance.   Thus,  systems  do  not  respond 
immediately  to  inputs,  and  this  must  be  considered  in 
the  system  design.   The  response  of  these  circuits  has 
the  following  appearance. 
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The  time  required  for  the  response  to  rise  to  63% 

of  the  maximum  is  known  as  the  time  constant  of  the  system 

The  narrower  the  resonance,  that  is,  the  narrower  the 

filter  bandpass,  the  longer  the  time  constant.   Thus, 

if  we  desire  to  resolve  frequencies  to  0.1  Hz,  we  must 

allow  enough  integration  time  for  the  filter  to  respond 

to  each  frequency,  and  this  will  be  longer  than  the  time 

required  to  resolve  frequencies  to  1  Hz. 
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Another  consideration  in  the  design  of  active  ranging 
detection  systems  is  the  requirement  to  accomodate  doppler 
in  the  returning  echoes.   The  input  filter  must  have  a 
bandwidth  wide  enough  to  pass  echoes  that  are  doppler 
shifted  by  the  maximum  amount  anticipated  during  operation 
of  the  system,  or  the  echoes  will  not  be  detected.   If 
a  filter  of  sufficient  bandwidth  to  ensure  this  does  not 
give  a  sufficient  S/N  ratio  improvement,  one  may  alter- 
natively use  a  bank  of  narrow  band  filters  centered  at 
various  frequencies  corresponding  to  different  doppler 
shifts.   This  is  illustrated  in  the  diagram  on  the  next 
page. 
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V.   RANDOM  SIGNALS,  POWER  SPECTRAL  DENSITY,  AND  NOISE 

A.   RANDOM  SIGNALS 

Heretofore  the  signals  discussed  have  been  explicitly 
described  "deterministic  signals",  such  as  a  square  wave, 
sinusoidal  wave,  sawtooth,  etc.   In  describing  a  signal  ex- 
plicitly as  sume  function  of  time,  it  is  assumed  that  its 
amplitude  and  phase  are  known  exactly  for  all  time;  i.e., 
past,  present,  and  future.   No  real  world  signal  meets  these 
criteria,  for  then  the  signal  would  convey  no  new  informa- 
tion, as  it  could  be  predicted  exactly  for  any  future  time. 
Therefore,  consideration  must  be  given  to  random  signals. 
Random  signals  of  interest  to  the  Antisubmarine  Warfare  spe- 
cialist emanate,  for  example,  from  various  machinery  and  flow 
noises.   Noise  in  the  ocean  due  to  waves,  wind,  biological 
sources,  seismic  disturbances,  etc.,  is  a  special  case  of 
random  signal,  and  it  is  always  present  to  some  degree  wheth- 
er  or  not  a  "target"  is  present.   Processing  equipment  must 
be  able  to  distinguish  between  the  random  signal  of  interest, 
unwanted  signals,  and  noise. 

A  random  signal,  unlike  a  deterministic  signal,  may  only 
be  described  by  its  statistical  characteristics.   Its  past 
history  and  its  predicted  future  will  be  stated  in  terms  of 
some  average  value  and  a  given  probabi lity  that  the  signal 
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will  be  within  certain  limits  at  a  specified  time.   This 
"given  probability"  is  represented  by  a  probability  density 
function  (pdf ) . 

What  is  a  probability  density  function?   A  function  is  a 
descriptor,  i.e.,  it  usually  describes  something.   Density  is 
a  measure  of  quantity  contained  within  a  specified  space. 
Probability  refers  to  the  chance,  or  likelihood,  of  the  occur- 
rence of  an  event.   Therefore,  a  probability  density  function 
describes  the  likelihood  of  some  quantity  being  contained 
within  a  specified  space. 

How  does  the  probability  density  function  relate  to  a 
random  signal?   Assume  a  random  signal,  x(t),  which  has  some 
randomly  varying  voltage  amplitude  over  time. 
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Plot  the  signal  with  a  vertical  scale-  representing  voltage 
amplitude   and  the  horizontal  scale  representing  time. 
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Suppose  now  that  l~volt  represents  any  value  between  0.5  and 
1.49;  2-volt  represents  J..5  to  2.49;  3-volt  represents  2.5  to 
3.49,  etc.   From  the  random  signal,  a  plot  is  now  made  with 
the  vertical  scale  still  representing  voltage  amplitude  but 
with  the  horizontal  scale  representing  the  number  of  times 
the  signal  is  at  that  amplitude,  or  within  that  amplitude 
range. 

During  time  interval  1,  the  signal  was  present  only  in 
the  3-volt  range  and  that  one  block  is  shaded  in  on  the  plot. 
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During  time  interval  2,  the  signal  is  present  in  the  1, 
2,  and  3-volt  ranges.   Those  additional  blocks  are  shown  in 
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the  next  diagram. 


1C^ 


D 


m 


1  '  2  '  3  '  4-  '  S  I    i  ■  7  '  t  '  1  '  ic 
fl  U  »l  A  e  >-        cf'      l*  m*>     thi 


ft.  m 


i  / ."  r  n.  </e        r  tf-  /i  <?  <?• 


During  time  interval  3,   the  signal  is  present  in  tne  ]_ 
and  2-volt  ranges.   Those  additional  blocks  are  shaded  in. 
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The  development  should  start  to  become  apparent.   Look 
now  at  time  intervals  4, 5, and  6.   The  signal  is  present  in 
the  3-volt  range  once,  the  2-volt  range  twice,  and  in  the  1- 
volt  range  three  times.   Shade  those  in. 
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To  finish  the  plot,  look  at  time  intervals  7,8,9  and  10. 
The  signal  appears  in  the  4-volt  range  four  times;  in  the  3- 
volt  range  three  times;  in  the  2-volt  range  four  times;  in 
the  1-volt  range  four  times;  and  the  0-volt  range  once. 
Shade  them  in. 
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The  plot  above  represents  the  amplitude  ranges  that  the 
signal  was  in  for  the  entire  time  it  was  observed.   The  total 
number  of  shaded  blocks  represents  the  number  of  times  the 
signal  is  within  the  various  amplitude  zones  during  the  var- 
ious time  intervals  for  the  entire  duration  of  the  signal: 
the  29  blocks  account  for  100%  of  signal  time.  (Count  them.') 
Each  shaded  block  is  regarded  as  a  sample  value.   If  the  num- 
ber of  sample  values  is  denoted  by  n,  then  n  =  29  in  this 
example.   The  percentage  of  time  the  signal  is  in  a 
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particular  voltage  range  can  be  approximated  by  dividing  the 
number  of  signal  sample  values  in  that  voltage  range  by  the 
total  number  of  sample  values,  n. 

0-volt  range:   1/29  =   3% 
1-volt  range:   9/29  =31% 
2-volt  range:   8/29  =  28% 
3-volt  range:   7/29  =  24% 
4-volt  range:   4/29  =  14% 
Total   100% 
The  percentage  of  time  the  signal  spends  at  other  ranges 
can  also  be  calculated.   Suppose,  for  example,  that  in  an  ASW 
system  there  is  a  detector  which  detects  only  signals  with 
amplitudes  of  3  volts  or  more  (the  detector  has  a  "detection 
threshold"  of  3  volts).   To  determine  the  percentage  of  time 
the  above  random  signal  is  detectable,  it  is  necessary  to  de- 
termine the  percentage  of  time  it  remains  in  the  range  of  3 
or  more  volts.   The  3-volt  range  covers  from  2.5  to  3.49 
volts.   Take  one-half  the  time  the  signal  is  in  the  3-volt 
range  and  add  it  to  the  entire  time  the  signal  is  in  the  4- 
volt  range  to  approximate  the  total  time  the  signal  is  at  or 
above  3  volts. 

(1/2) (7)  +  4  =  7.5   Divide  by  n  to  get  percentage. 
7.5/29   =   26%     Therefore,  the  signal  is  detectable 

26%  of  the  time. 
The  accuracy  in  the  determination  of  the  percentage  of 
time  the  signal  is  within  any  specified  frequency  range  can 
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be  increased  by  making  the  voltage  amplitude  ranges  smaller 
and  shortening  the  duration  of  the  time  intervals.   The 
next  diagram  shows  sample  values  for  the  same  signal,  but 
with  amplitude  ranges  of  1/2-volt  and  the  time  interval  dura- 
tion halved.   It  is  constructed  in  the  same  amnner  as  the 
previous  plot.   A  check  to  verify  its  accuracy  will  help  you 
evaluate  your  understanding  of  the  procedure. 
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As  amplitude  range ' and  the  time  interval  length  are  made 
smaller  and  smaller,  the  plot  approaches  a  continuous  curve. 


continuous  curve 
approximation 
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In  most  practical  applications  a  formula  for  the  continuous 
curve  can  be  determined,  and  the  corresponding  function  inte- 
grated over  the  amplitude  range  of  interest  to  find  the  per- 
centage of  time  the  signal  was  within  that  range. 

The  probability  density  function  of  a  signal  represents 
the  "probability"  that  a  signal  will  be  within  an  amplitude 
range.   As  the  plot  constructed  above  represents  100%  of  the 
time  of  the  signal,  so  too  does  the  probability  density  func- 
tion represent  the  total  probability  of  the  signal.   It  is 
scaled  from  0  to  1,  i.e.,  the  area  under  the  probability  den- 
sity function  curve  is  equal  to  ]_.       If  the  detector  in  the 
example  just  mentioned  detects,  with  certainty,  any  voltage 
amplitude  of  3  volts  or  more,  and  the  signal  spends  26%  of 
the  time  at  3  volts  or  more,  then  the  probability  that  the 
signal  will  be  detected  in  the  3-volt  or  more  range  is  0.26. 
The  probability  density  is  expressed  by  a  function  which  may 
be  integrated  to  find  the  probability  that  the  signal  is  in 
any  prescribed  amplitude  range.   For  probability  density  func- 
tions used  in  many  applied  problems,  the  integration  process 
can  be  avoided,  sincd  the  results  are  readily  available  in 
standard  tables. 

The  probability  density  function  (pdf)  most  commonly  used 

for  electrical  signals  is  the  "gaussian"  or  "normal"  pdf.   It 

is  the  familiar  bell-shaped  curve  described  by  the  equation 

-(x-m)2 


P  (x)  = 


1    ~  2a2 


/2tt 
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In  the  equation  for  p (x) , 

m  =  mean,  or  average,  value, 

i 

x  =  value  of  the  signal,  x(t),  at  time  t, 

a  =  standard  deviation  (o  is  the  Greek  lower  case  sigma) . 
The  standard  deviation,  a,  is  a  measure  of  the  spread  of  the 
signal  abour  its  average  value,  m.   In  the  next  figure  we 
we  show  two  gaussian  curves  with  the  same  mean,  but  different 
values  of  a.   The  area  under  each  curve  is  equal  to  1.0. 


wo 

For  the  gaussian  pdf,  about  6  8%  of  the  area  under  the  curve 
lies  between  m  -  a  and  m  +  a,  which  means  that  the  probability 
that  the  random  signal,  x(t),  assumes  a  value  between  m  -la 
and  m  +  \o    is  0.68. 


m 


+ 


la,  P(x)  =  0.6: 


(Shaded  area  =6  8% 
of  total  area) 

m  -    2a,    P(x)    =    0.95 
(Shaded   area   =    95% 
of   total    area) 


:'Ler       r-n  +2<7~ 


81*. 


m  -    3a,    P(x)    =    0.997 
(Shaded   area  =    99.  7% 
of   total    area) 


-3  0~ 
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In  the  illustrations  on  the  previous  page,  only  0.99  7  of 
the  total  probability  has  been  accounted  for.   This  is  because 
there  are  no  limits  on  a  gaussian  pdf;  the  range  is  unbounded. 
A  signal  may  conceivably  exist  at  any  finite  amplitude  above 
or  below  the  average  value.   The  probability  that  the  signal 
will  occur  outside  the  m  -  3a  range  is  only  0.003  (1.0  -  0.997 
=  0.00  3).   Limiting  circuits  are  normally  used  early  in  the 
processing  stage  to  prevent  the  rare  occurrence  of  a  large 
amplitude  signal  from  damaging  equipment.   The  important  part 
of  the  signal,  0.997  of  it,  is  still  retained. 

The  gaussian  pdf  is  centered  around  an  average  value. 
For  an  electrical  signal  this  is  the  average  value  of  the 
amplitude  of  the  signal.   Look  at  the  signal  and  its  pdf 
side  by  side  with  the  pdf  oriented  to  the  signal  amplitude. 
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The  value  chosen  for  m  is  the  statistical  average  of  the 
29  sample  values  above.   Each  of  the  sample  values  was  0,  1, 
2,  3,  or  4.   The  average,  m,  is  equal  to  the  sum  of  the  indi- 
vidual sample  values  divided  by  the  sample  size,  n. 

0+1+1+1+1+1+1+1+1+1+2+2+2+2+2+2+2+2+3+. . .+4. . . 

29 

or,  since  many  values  are  the  same, 

1(0)  +9  (1)  +8(2)  +  7(3)  +  4(4)       , 

29  l 

-,   n  n 

Algebraically  expressed,  this  is  m  =  —  .1,  x.  where  the  .£, 

means  the  sum  of  all  the  samples,  x.,  from  j  =  1  to  j  =  n,  or 

x,  +  x0  +  x-.  +.  .  .  +  x  .   If  the  number  of  sample  values  becomes 
1     2     3        n  c 

very  large,  the  sum  of  the  products  of  the  individual  sample 
values  and  their  probability  of  occurrence  becomes  a  better 
choice  for  m.   Increasing  the  number  of  sample  values  gives  a 
more  accurate  description  of  the  signal.   If  m  is  calculated 
for  the  plot  with  the  amplitude  range  and  time  duration  halved, 
n  =  74   and  m  =  2.095.   The  algebraic  expression,  if  n  is  very 

large,  is  approximately  °° 

m  =  .Z,  (x. )•  P (x. ) 
D=l    3  3     • 

It  has  been  shown  how  statistical  characteristics  may  de- 
scribe a  signal.   The  next  problem  is  how  to  obtain  signal 
characteristics  that  are  statistically  valid  to  identify  a 
signal. 

The  validity  of  statistical  methods  in  determining  the 
outcome  of  a  process  depends  on  repition  of  the  process  over 
and  over  again.   With  a  random  signal,  x(t),  this  repetition 
may  be  accomplished  by  using  an  "ensemble"  of  random  signal 

86 


generators  with  each  signal  generated  having  the  same  broad 
statistical  characteristics.   If  the  individual  outputs  from 
this  "ensemble"  are  measured  simultaneously  at  some  time,  t, 
and  averaged,  the  result  would  be  an  "ensemble  average",  in- 


dicated by   x(t).   This  "ensemble  average"  is  identical  to, 
and  has  all  the  properties  of,  a  statistical  average.   How- 
ever, most  signals  of  interest  in  an  ASW  environment  originate 
from  a  single  source.   Therefore,  only  the  time  average  of 
the  signal,   ^x(t)\  ,  may  be  obtained. 

The  time  average,  ^(tn,  and  the  ensemble  average,  x(t), 
are  not  always  the  same.   For  example,  suppose  the  statisti- 
cal characteristics  of  the  signal  from  the  random  signal  gen- 
erators in  the  ensemble  vary  with  time.   Such  a  variation 
would  not  be  reflected  in  measurements  made  at  a  fixed  time, 
t,.   The  ensemble  average  at  time  t-.  would  be  different  from 
the  ensemble  average  at  time  t^ .   A  random  signal  with  these 
characteristics  is  called  "non-stationary".   If  a  detector  is 
used  to  detect  the  random  signal  in  our  example  above,  what 
is  the  effect  on  the  detector  if  the  random  signal  is  "non- 
stationary"?   The  average,  m,  was  calculated  as  2.14.   Suppose 
that  m  varies  with  time,  and  that  over  another  time  interval 
m  is  calculated  to  be  L. 14.   If  the  signal  over  the  second 
time  has  the  same  characteristics  as  over  the  first,  the 
effect  on  detection  would  be  the  same  as  raising  the  thresh- 
old to  4  volts  over  the  first  interval,  as  indicated  in  the 
next  diagram. 
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The  probability  that  x  =  3  volts  or  more  would  become  2/29  = 
0.069  vice  0.26  for  m  =  2.141   If  the  average  is  unchanged 
and  the  standard  deviation  increased  by  a  factor  of  three, 
the  detector  would  be  observing  only  0.68  of  the  signal  as 
opposed  to  0.99  7.   Therefore,  if  the  random  signal  is  not 
"stationary"  (i.e.,  its  statistical  characteristics  do  not 
vary  with  time)  or  at  least  "stationary"  for  the  processing 
time  of  the  detector,  the  detector  will  be  very  inefficient. 
Even  though  the  signal  may  be  stationary,  the  time  and  ensem- 
ble averages  may  still  be  unequal. 

Fortunately,  many  random  signals  that  mjst  be  processed 
in  ASW  have  time  and  ensemble  averages  which  are  identical, 
perhaps  not  for  all  time,  but  at  least  for  the  processing 
time  required.   Signals  with  their  time  and  ensemble  averages 
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equal,  ^x(t)^  =  x(t),  are  said  to  come  from  an  "ergodic  pro- 
cess",or  to  have  the  property  of  "ergodicity" .   An  "ergodic 
process"  is  also  stationary.   Under  the  assumption  that  a 
signal  from  an  ergodic  process  means  the  time  and  ensemble 
averages  are  equal,  the  signal  is  stationary  and  the  time 
average  will  have  all  the  properties  of  a  statistical  average 

Having  extablished  that  a  time  average  will  be  statisti- 
cally valid  if  it  is  from  an  ergodic  process,  how  do  we  ob- 
tain the  other  statistical  characteristics  which  describe  a 
signal?   All  acoustic  detectors  convert  the  signal  (and,  un- 
fortunately, the  noise,  which  will  be  addressed  later)  into 
some  time- varying  voltage  or  current,  x(t). 


SCCfia- 
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The  time  average  of  the  signal,  <^(t)\  ,  assuming  ergodi- 
city, is  equal  to  the  average  value,  m.   This  average  value 
is  both  the  value  about  which  the  signal's  gaussian  pdf  is 
centered  and,  as  the  time  average,  represents  the  "dc"  compo- 
nent of  the  signal  (direct  current,  also  called  the  steady- 
state  or  the  non-time-varying  component).   If  m  =  0,  there 
is  no  "dc"  component.   If  the  signal,  x(t),  is  impressed  a- 
cross  a  one-ohm  resistor,  the  instantaneous  power  dissipated 
is  x2(t).   The  time  average,  /;2(t)\,  becomes  the  "second 
moment",  which  represents  the  total  average  power,  x2 ,  of  the 
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electrical  signal.   The  "second  moment"  differs  from  the 
"average  value  squared",  m2,  which  is  the  dc  power,  or  the 
power  in  the  non-time-varying  component:  m2  is  obtained  by 
squaring  after  averaging,  whereas  x2  is  obtained  by  squaring 
before  averaging. 

As  an  illustraion  of  these  values,  consider  a  random 


signal,  x (t) . 
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The  following  drawings  illustrate  computation  of  m2  and 


x2. 
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There  sre  two  other  measures  which  identify  a  specific 
gaussian  pdf  besides  the  average  value,  m,  about  which  it  is 
centered.   A  random  signal  is  composed  of  a  dc  component  and 
an  "ac"  component  (alternating  current,  or  the  time  varying 
component  of  the  signal).   The  total  power,  x2 ,  must  be  the 
sum  of  the  dc  power,  m2 ,  and  the  ac  power,  a2.   Though  repre- 
senting the  ac  power  of  the  signal,  a2  is  the  "variance"  of 
the  gaussian  pdf.   The  variance  is  equal  to  the  second  moment 
minus  the  average  squared,  a2  =  x2  -  rn2.   For  an  electrical 
signal  the  variance  is  the  total  average  power  minus  the  dc 
power.   As  the  standard  deviation,  a,  represents  the  square 
root  of  the  ac  power,  it  is  more  commonly  known  as  the  root- 
mean  square  (rms)  value.   We  summarize  these  definitions 
briefly. 

Given  an  ergodic  random  signal,  x(t): 

1.   the  average  (mean)  value,  m,  is  its  dc  component; 
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2.  the  average-squared,  m2 ,  is  its  dc  power; 

3.  the  second  moment,  x2  ,  is  its  total  average  power; 

4.  the  variance,  a2  =  x2  -  m2 ,  is  its  ac  power; 

5.  the  standard  deviation,  o,    is  its  rms  value. 

In  "Fourier  Transform  Properties",  the  concept  of  correla- 
tion of  a  signal  was  discussed.   It  is  considered  again  here 
because  of  its  unique  properties,  and  because  of  its  relation- 
ship to  the  "power  spectral  density"  of  a  random  signal. 
Those  who  feel  they  have  not  grasped  the  concept  of  correla- 
tion should  review  that  section  before  proceeding* 

The  autocorrelation  function,  R(t),  was  defined  previously 
as  the  time  average,  R(t)  =  ^x(t).-x(t  +  t)>.   For  an  ergodic 
random  signal  the  autocorrelation  function  becomes: 

R(t)  =  ^t(t)  .  x(t  +  t)\   =  x(t)  .  x(t  +  t) 
The  value  of  R(t)  is  a  function  of  the  time  delay,  t.   At  t=0  , 
R(t)  is  a  manimum: 


R(0)  =  x(t)»  x(t+0)  =  x  =  total  average  power. 

As    increases,  R(t)  decreases,  until  at   x=  infinity,  R(t) 
is  a  minimum: 

R(»)  =  m2  =  dc  power; 

R(0)  -  R(°°)  =  a2  =  x2  -  m2  =  ac  power. 
The  autocorrelation  function  is  also  a  measure  of  the 
"time  coherence"  of  the  random  signal.   If  t  is  small,  x(t) 
and  x(t+x)  will  be  close  together  in  time.   Therefore,  the 
presence  of  one  will  have  some  effect  on  the  presence  of  the 
other,  a  condition  known  as  "statistical  dependence".   If  the 
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presence  of  x(t)  has  no  effect  on  x(t+i),  such  as  happens  if 
t  is  large  and  x(t)  is  random  and  non-periodic7  then  x(t)  and 
x(t+x)  are  said  to  be  "statistically  independent".   If  for  a 
signal  with  a  gaussian  pdf  the  absolute  value  of  the  auto- 
correlation function  minus  the  dc  power  equals  zero,  |R(t)|  - 
m2  =  0 ,  then  x(t)  and  x(t+x)  are  statistically  independent 
and  "uncorrelated" .   Statistical  independence  simplifies  the 
calculation  of  the  autocorrelation  function.   Whether  or  not 
two  signals  are  "correlated"  becomes  important  in  the  calcula- 
tion of  their  average  powers. 

Assume  that  a  signal  has  been  formed  by  the  addition  of 
two  separate  signals:  z(t)  =  x(t)  +  y(t).  The  correlation 
function  of  z (t)  will  be  of  the  form 

(x+y) (x+y)  =  x2  +  xy  +  yx  +  y2 ,   and  is  equal  to: 


R 
z 


(t)  =  <x(t)»  x(t+T^  +^(t)-  y(t+x)^  +  ^(t).y(t+x))  +  ^  (t)  •  x  (t+T  )> 


Rz(t)  -    Rxx(T)     +   Ryy(T)     +   Rxy(T)     +   RYx(t) 
where  the  last  two  terms  are  the  "cross  correlation  terms"  of 

x(t)  with  y(t)  and  y(t)  with  x(t).   As  noted  previously,  the 

above  correlation  functions  represent  the  average  power  in 

the  signals,  and  can  be  written  as: 

p=p     +P     +   P     +P. 
z    xx      yy      xy      yx 

If  the  two  signals,  x(t)  and  y(t),  are  correlated  over  the 
time  delay,  T,    then  the  total  signal  power,  S,  is  the  sum  of 
the  signal  power  in  signal  x(t),  Sxx,  the  signal  power  in  sig- 
nal v(t),  S   ,  and  the  cross  correlation  of  their  combined 

J     yy 

signal  powers,  S    and  S   : 

^.y      yx 
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s.=  s   +s   +  s   +  s 

z    xx    yy    xy    yx 

If  the  signal  is  the  sum  of  mutually  un correlated  signals, 

as  in  the  case  of  noise  uncorrelated  over  the  time  delay,  t, 

then  the  total  noise  power,  N,  is  the  sum  of  only  the  noise 

power  N    and  N   .   The  cross  correlations  are  equal  to  zero: 
j^j^  yy  '       ■  ■  "  ' 

N   =  N    +  N   ;   (N    =0;   N    =  0)  . 
z     xx     YY      xy         yx 

For  a  signal  which  is  the  sum  of  mutually  uncorrelated 
signals,  the  correlation  function  is  the  sum  of  only  the  auto- 
correlations (the  cross  correlations  equal  zero) ,  and  the  av- 
erage power  is  the  sum  of  only  the  average  powers  of  the 
individual  signals  (the  cross  correlated  powers  equal  zero) . 

B.   POWER  SPECTRAL  DENSITY 

j.iius  .Lai.  on  xy  iflc  uimc  QGiTiaj-ri  3sp£Ci.  Oi  random  signals  uSG 
been  addressed.   The  frequency  domain  is  represented  by  the 
power  spectral  density,  which  was  introduced  in  "Filters  and 
Bandwidth  Requirements".   The  device  used  to  transfer  between 
the  time  and  frequency . domains  is  the  Fourier  Transform.   The 
Weiner-Kinchine  Theorem  states  that  the  power  spectral  density, 
G(f),  and  the  autocorrelation  function,  R(x),  are  Fourier 
Transforms  of  each  other: 

G(f)  =  J    [R(t)]   =   fx(t)  x(t+T)e_:)a3TdT  . 

And,  from  the  Duality  Theorem, 

if  R(t)  «~K3(f),   then   G(t)<-+  R(f). 

As  can  be  seen  above,  the  random  signal  may  be  described 

in  either  the  time  domain,  R(t),  or  the  frequency  domain,  G(f), 
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with  free  interchange  between  the  two. 

Since  the  power  spectral  density  is  the  Fourier  Transform 
of  the  autocorrelation  function,  the  Fourier  Transform  Theo- 
rems may  be  applied  to  the  power  spectral  density.   The  sys- 
tems designer  uses  the  properties  expressed  by  these  theorems 
to  design  an  efficient  detector/signal  processing  system.   The 
ASW  specialist  who  has  knowledge  of  these  properties  can  bet- 
ter understand  how  his  system  is  designed  to  operate  and 
therefore,  utilize  it  more  effectively.   The  duality  theorem 
was  already  utilized.   The  other  major  theorems  applicable  to 
power  spectral  density  will  be  explained  along  with  examples 
of  their  usage. 

]_.   Frequency  Translation  Theorem 

If  a  signal,  x(t),  is  bandlimited  in  W  <  f  ,  and  a 
new  signal,  y(t),  is  formed  by  multiplying  x(t)  by 
cos  (a)  t+<J>)  ,   y(t)  =  x(t)  cos  ( w  t+cj>)  ,  then  the  correlation 
function  of  y(t)  is  equal  to  one-half  the  correlation 

function  of  x(t)  cosco  t,  that  is  R  (t)  =  (]./2)R  (t)cosco  t. 

c  y  x       v- 

In  addition,  the  power  spectral  density,  G(f),  is 
Gy(f)  =  (l/4)Gx(f-fc)  +  (l/4)Gx(f+fc)  • 

Multiplication  in  the  time  domain  by  cosco  t  becomes  transla- 
tion in  the  frequency  domain,  as  the  power  spectrum  is  shifted 
up  and  down  in  frequency  by  ±f  .   The  condition  that  the  band- 
width, W,  be  less  than  f   is  necessary  to  ensure  that  Gx(f-fc) 

and  G  (f+f  )  do  not  overlap,  otherwise  aliasing,  as  described 
x    c 

in  "Filters  and  Bandwidth  Requirements"  would  result.   Notice 
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also  that  in  multiplication  by  cos  (oj-  t+4>)  ,  the  phase  factor, 
<J>,  is  lost  in  translation  to  the  power  spectral  density. 
Therefore,  the  original  random  waveform  may  not  be  recon- 
structed from  knowledge  of  the  power  spectral  density.   How- 
ever, by  returning  to  the  time  domain  via  the  inverse  Fourier 
Transform,  the  autocorrelation  function,  R(t),  with  its  sig- 
nificant statistical  properties,  may  be  found. 

Frequency  translation  is  very  useful  in  signal  processing 
Suppose  that  due  to  space,  weight,  or  power  limitations  it  is 
impractical  to  locate  the  detector  and  the  processor  together, 
The  low  frequency  (0  to  ]000  Hz)  sound  may  be  detected,  trans- 
lated to  radio  frequency  (RF)  (HF  =2-30  MHz,  VHF  =  100  - 
200  MHz,  UHF  =  200  -  400  MHz)  and  broadcast  to  the  processor. 
In  the  same  manner,  one  processor  can  handle  many  detectors 
if  each  detector's  signal  is  translated  to  a  separate   fre- 
quency band  within  the  RF  range  (Detector  ].,  frequency  f-,  to 
f ? ;  Detector  2,  frequency  f~  to  f .;    etc.,  all  within  the  UHF 
band) .   The  problem  could  be  one  of  processing  signals  of 
many  different  frequencies.   Rather  than  having  many  differ- 
ent processors,  each  processing  only  one  signal,  it  may  be 
better  to  have  one  processor  operating  at  a  fixed  frequency. 
All  of  the  signals  could  be  translated  to  the  processor  fre- 
quency for  processing.   If  the  processor  is  very  elaborate 
and/or  very  expensive,  this  method  becomes  desirable. 

2.   Integration  Theorem 

rr  i 

If  y(t)  =/  x(t')dt\   then  G  (f)  =  ± Gv(f)- 

<co  y      (27Tf)2    x 
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The  multiplication  factor,  ,_  _. 2  ,    indicates  that  the  low 

(27Tf)  

frequency  components  of  the  power  spectral  density  will  be 
enhanced  by  integration.   The  ocean  acts  as  an  "integrating 
filter",  in  that  it  allows  the  low  frequencies  to  propagate, 
while  rapidly  attenuating  the  higher  frequencies.   This  is  one 
reason  why,  in  an  attempt  to  get  long  detection  ranges,  active 
sonars  and  passive  processing  systems  have  increasingly  ex- 
ploited the  lower  frequencies. 
3.   Differentiation  Theorem 

This  theorem  is  offered  without  an  example  since  it 
follows  so  naturally  the  integration  theorem. 

If  Y(t)  =  d[x(t)]/dt,  then  G  (f)  =  (2^f ) 2G  (f). 

y  x 

The  multiplication  factor,  (2-rrf)2,  indicates  that  the  high 
frequency  components  of  the  power  spectral  density  will  be 
enhanced  by  differentiation. 

C.   NOISE 

At  the  beginning  of  this  section,  noise  was  mentioned  as 
a  "special  case  of  random  signal  which  is  always  present  to 
some  degree  whether  or  not  a  'target'  is  present".   In  com- 
munications, noise  is  often  defined  as  any  electrical  inter- 
ference, or  unwanted  signal.   Signal  then,  in  addition  to  the 
characteristics  previously  used  to  describe  it,  has  the  qual- 
ity of  being  wanted,  sought  after,  emanating  from  a  source 
of  interest,  and  containing  intelligence.   In  ASW,  signal  is 
always  associated  with  "target".   Noise,  on  the  other  hand, 
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is  anything  which  interferes  with,  or  tends  to  mask,  the 
signal.   Often  "one  man's  signal  is  another  man's  noise"  as 
any  patrol  aviation  Tactical  Coordinator  whose  buoy  field  has 
just  been  penetrated  by  a  "pinging"  Destroyer  will  verify. 
An  acoustic  detector  converts  both  signal  and  noise  into  a 
time-varying  electrical  voltage  or  current.   An  examination 
of  some  of  the  characteristice  of  noise  will  help  in  under- 
standing how  systems  may  be  designed  to  reduce  it. 

Another  way  to  define  noise  is  simply  as  that  quantity 
observed  in  the  absence  of  signal.   Noise  may  come  from  a  var- 
iety of  manmade  or  naturally  occurring  sources.   Systems  are 
designed  to  eliminate  as  much  noise  as  possible,  but  some 
noise  will  inevitably  remain.   The  thermal  motion  of  electrons 
in  the  conducting  media  ot  the  components  of  a  system,  for 
instance,  is  one  unavoidable  cause  of  electron  noise.   Ther- 
mal noise  is  interesting  in  that  its  power  spectrum  is  con- 
stant over  a  wide  range  of  frequencies.   It  is  designated 
■**white  noise"  by  analogy  to  white  light,  as  all  frequency  com- 
ponents are  present  in  equal  amplitudes.   The  amplitude  of 
thermal  noise  has  been  proven  to  have  a  "gaussian"  pdf.   It  is 
therefore  referred  to  as  "gaussian  white  noise"  ("gaussian" 
amplitude  distribution;  "white"  frequency  distribution) .   If 
"gaussian  white  noise"  is  played  over  a  loudspeaker,  it  sounds 
dull  and  monotonous,  somewhat  like  a  waterfall.   The  subtle- 
ties of  its  random  variations  are  hidden  from  the  human  ob- 
server. 

The  power  spectrum  of  gaussian  white  noise,  with  zero 
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mean,  is  G(f.)  =  —  ,  where  n  is  power  density  in  watts  per 
Hertz.   By  Fourier  transformation,  the  autocorrelation  func- 


tion   is 


R(t) 


JUT 


df  =    (^)6(t)    ,    whe: 


/«\w     ,    wnere    6    represents 
the    "delta    function"    introduced  previously. 
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It  is  apparent  from  the  R(t)  plot  that  for  any  time  delay,  t, 
other  than  zero,  the  autocorrelation  function  is  zero;  there- 
fore, any  two  different  samples  of  a  gaussian  white  noise  sig- 
nal are  uncorrelated  and  statistically  independent. 

The  signals  of  interest  in  ASW  are  frequency  limited, 
i.e.,  there  is  some  finite  frequency  band  within  which  the 
signal  exists.   Filtering  in  order  to  improve  detection  of 
that  signal  was  discussed  earlier.   If  gaussian  white  noise 
is  filtered,  its  frequency  components  are  naturally  those  of 
the  filter  bandwidth,  but  the  amplitude  distribution  remasins 
gaussian.   The  output  power  spectrum  of  white  gaussian  noise, 
filtered  by  an  ideal  filter,  is  a  rectangular  function, 

X]  f 

Gy^f)  ^~2"'  ^2B)  >  where  II  indicates  a  rec- 
tangular function  and  B  represents  the  filter  bandpass.  The 
autocorrelation  function  is  the  inverse  Fourier  Transform  of 
the  power  spectral  density.  As  indicated  in  "Fourier  Trans- 
form Properties",  the  Fourier  Transform  of  a  rectangular 
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The  above  figure  shows  that  filtering  has  produced  the  follow- 
ing results. 

1.  The  power  spectrum,  though  no  longer  white,  is  con- 
stant over  the  finite  frequency  range  of  the  filter. 

2.  The  output  power  is  finite,  y2  =  j    (2B)  =  nB.   It 
varies  linearly  with  bandwidth,  B. 

3.  The  output  signal  is  correlated  over  time  intervals 
1 


of  about 


2B  * 


How  is  this  knowledge  used  in  signal  processing? 

Item  (1).   The  noise  in  which  a  signal  of  interest  must 
be  detected  is  often  not  gaussian  white  noise.   However, 
through  judicious  choice  of  a  bandwidth,  its  power  spectral 
density  may  be  constant  over  that  range.   Otherwise,  a  pre- 
whitening  filter  may  be  used  to  make  the  noise  power  spectral 
density  constant  over  the  bandwidth.   Therefore,  after  filter- 
ing it  may  be  treated  as  gaussian  white  noise. 
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Item  (2).   The  noise  power  after  filtering  is  y2  =   nB. 
It  is  a  direct,  linear  function  of  bandwidth,  B.   By  narrow- 
ing the  bandwidth,  the  amount  of  noise  power  present  is  re- 
duced.  Therefore,  the  signal  power  is  enhanced  in  relation  to 
the  noise  power.   This  is  crucial  in  a  power,  or  energy 
detector. 

Item  (3).   The  bandwidth  of  the  filter  determines  the  min- 
imum time  delay,  t ,  in  a  correlation  detector.   It  must  be 
greater  than  ■—,     (t>  -k^)  ,    or  the  noise,  as  well  as  the  signal, 

Z.  B         Z.  B 

will  be  correlated.   If  the  noise  is  uncorrelated,  which  it 
will  be  for  t>  ■—,    only  its  average  power  will  add  in  a  corre- 

Z.  B 

lation  detector,  thus  decreasing  in  relation  to  the  signal 

power. 

The  bandwidth  of  the  filter  also  determines  the  finite 

integration  time  to  be  used  in  determining  the  time  average. 

The  time  average  is  indicated  by  \j (brackets)  and  is  defined 

as  t  i   >»i 

lim  £  /  x(t)  x(t+x)  dt. 

Replacing  T  by   ^      and     by   —^,    we    obtain 

1_ 

bJ2B      x(t)    x(t+I-)dt. 

A  review  of  the  development  of  bandwidth  versus  integration 
time  in  the  "Filters  and  Bandwidth  Requirements"  section  may 
be  helpful  in  the  understanding  of  these  relationships. 

This  section  has  been  an  attempt  to  describe  the  concepts 
of  random  signals,  power  spectral  density,  and  noise.   The 
intent  has  been  to  provide  a  basis  for  understanding  their 

101 


individual  characteristics.   A  grasp' of  these  characteris- 
tics is  necessary  to  understand  the  methods  of  signal  process- 
ing.  The  most  common  methods  in  current  use,  DELTIC,  energy 
detectors,  correlation  detectors  and  beam  formers,  will  be 
presented  in  following  sections. 
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VI.  THE  DELTIC 


In  the  section  on  filters  and  bandwidth,  the  relation- 
ship between  integration  time  and  frequency  resolution 
was  discussed.   It  was  shown  that  one  must  integrate  for 
longer  periods  of  time  as  the  frequency  resolution 
requirements  become  stricter.   In  many  cases,  however, 
the  integration  time  required  to  give  the  necessary 
resolution  is  much  too  long  to  be  of  tactical  use.   In 
order  to  circumvent  this,  several  processing  schemes 
have  been  developed  which  in  effect  give  more  resolution 
for  a  given  integration  time.   One  of  the  more  common 
in  use  with  military  processors  is  known  as  the  DELTIC 
(DELay  Time  Compressor) .   Let  us  investigate  the  manner 
in  which  this  processor  accomplishes  this  feat. 

It  has  been  shown  that  processing  the  signal  as  it 
is  received  by  the  detector  will  not  give  any  greater 
resolution.   The  DELTIC  takes  the  signal  as  received, 
and  records  it  on  a  magnetic  tape,  disc,  or  drum.   It 
does  this  in  a  special  sequence,  however.   To  see  this, 
let  us  take  a  simple  signal  and  observe  how  it  is  processed 


t-** 
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For  simplicity  in  showing  the  operation,  an  8-sample 
point  "time  window"  will  be  used.   In  actual  practice, 
1000  or  more  sample  points  are  commonly  used.   If  the 
time  between  sample  points  is  T,  the  time  window  observed 
will  be  8t   =  T.   Therefore  if  the  DELTIC  is  not  used, 
the  best  resolution  we  can  get  is  1/8T,  and  it  would 
require  8t  seconds  to  get  this  resolution.   The  physical 
layout  of  a  DELTIC  is  shown  in  the  next  diagram. 

INPUT  R£M>*C)UT 


The  drum  is  rotated  at  a  rate  such  that  it  makes  one 
revolution  plus  one  sample  space  during  the  time  period  t 
In  this  case,  with  8  sample  points  it  makes  1.125  rev- 
olutions per  t   seconds.   The  read-out  head  continuously 
reads  out  the  sample  values  as  the  drum  rotates  under  it. 
The  input  gate  injects  a  new  sample  point  once  each 
1.125  revolutions  in  such  a  way  that  the  sample  point 
which  has  been  on  the  drum  the  longest  time  is  replaced. 
Starting  with  a  blank  drum,  let  us  examine  the  process 
as  it  occurs. 
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When  the  drum  has  made  1.125  revolutions  (in   sec.) 
the  input  gate  inserts  sample  point  #2. 
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After  1.125  more  revolutions,  the  3rd  sample  point 
is  inserted. 


ZH^zirJjJ_Ll_Ll 


105 


The  remaining  spots  are  filled  similarly  until  all 
8  are  filled.  The  output  during  one  revolution  of  the 
drum  then  is: 


^Ajdnuum] 


When  the  next  sample  point  is  inserted  after  1.125  revo 
lutions,  it  displaces  the  first. 
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Once  the  drum  is  filled  (which  takes  8T  sec),  a  new 
value  replaces  the  oldest  one  every  T  seconds.   If  we 
look  at  the  output  from  the  read-out  head,  we  see  that 
there  are  9  values  read  out  each  T  seconds  (8  values/ 
revolution,  1.125  revolutions/r  seconds). 
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The  output  of  the  DELTIC  each  T  seconds  is  thus  the 
same  as  the  inputs  over  the  last  8T  seconds.   The  DELTIC 
operation  described  above  can  be  illustrated  by  the 
following  analogy.   If  we  record  80  seconds  of  conversation 
on  tape,  and  then  speed  up  the  playback  so  that  the 
conversation  takes  only  10  seconds,  we  have  compressed 
the  conversation  just  as  a  DELTIC  would.   The  same 
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information  is  in  the  10  second  playback  as  was  in  the 
80  second  conversation,  but  it  only  takes  l/8th  the  time 
to  hear  it.   We  also  notice  that  this  compressed  conver- 
sation has  been  shifted  in  frequency  by  the  compression, 
and  normal  voices  sound  like  Donald  Duck.   This  illustrates 
the  practical  effects  of  the  Fourier  scale  change  theorem. 

By  using  the  DELTIC  we  then  are  able  to  observe  the 
preceding  8T  seconds  worth  of  signal  during  each  interval 
of  length  f  seconds.   The  frequency  domain  plots  of  the 
original  signal  and  the  output  of  the  DELTIC  show  that 
the  original  spectrum  has  been  shifted  in  frequency  and 
stretched  out. 
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The  Fourier  scale  change  theorem  shows  that  we  have 
multiplied  the  time  domain  characteristics  of  the  signal 
by  the  factor  1/9,  and  so  the  frequency  domain  character- 
istics of  the  signal  have  been  divided  by  this  factor, 
thus  effectively  multiplying  the  frequency  plot  by  9. 

The  "compression  factor"  is  the  ratio  of  the  original 
sample  time,  T,  to  the  compressed  sample  time,  T'. 
Thus  in  this  case  it  is  the  ratio  T  to  T/9,  or  9.   Note 
that  not  only  is  the  center  frequency  and  bandwidth 
multiplied  by  this  factor,  but  also  the  minimum  resolvable 
frequency  difference.   Without  compression,  we  would 
integrate  for  T  seconds  to  obtain  a  resolvable  difference 
of  1/T.   Upon  compression,  this  becomes  9/T.   At  first 
glance,  this  does  not  seem  to  help,  since  the  resolution 
has  not  been  improved.   Why  use  the  DELTIC  if  it  doesn't 
help? 

The  benefit  in  using  the  DELTIC  is  that  we  can  obtain 
this  resolution  not  by  integrating  over  T  seconds,  but 
over  T/9  seconds  once  the  DELTIC  has  been  filled.   Many 
spectrum  analyzers  "scan"  the  frequency  band  of  interest 
by  using  a  single  band-pass  filter  of  bandwidth  equal 
to  the  minimum  resolvable  frequency  difference  of  the 
processor,  and  stepping  its  center  frequency  across  the 
band  one  step  each  integration  period.   Thus,  if  the 
filter  bandwidth  is  1/9  the  bandwidth  of  interest,  without 
the  DELTIC  it  would  require  9  integration  periods  (T)  for 
one  scan  of  the  band.   With  the  DELTIC,  each  integration 

108 


period  is  T/9 ,  and  thus  we  can  scan  the  entire  band  with 
the  same  resolution  in  what  was  previously  one  integration 
time,  T.   This  is  quite  an  advantage  for  a  band  of 
several  hundred  hertz  and  a  required  resolution  of  0.5  Hz. 

The  DELTIC  need  not  be  designed  for  insertion  of  a 
new  sample  value  every  revolution,  but  only  a  new  value 
approximately  every  10  revolutions.   In  this  way  one  can 
compress  the  signal  even  further  by  having  the  drum 
revolve  more  than  once  every  t   seconds. 

Instead  of  a  magnetic  drum,  disc,  or  tape  memory, 
a  digital  shift  register  may  be  used.   The  shift  register 
is  stepped  past  the  read-out  gate  in  such  a  way  as  to 
read  the  entire  register  each  r  seconds.   The  values 
are  then  recirculated.   Each  time  a  complete  circuit  is 
made,  a  new  value  is  input  to  replace  the  oldest,  and 
the  process  continues. 
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Some  systems  which  use  the  DELTIC  in  processing  are: 
a..   AN/AQA-7 

b.  AN/BQR-13 

c.  AN/SQS-26 

d.  AN/BQR-20 

The  illustrations  on  the  following  pages  show  a  side  by 
side  comparison  of  an  input  signal  and  the  resultant  DF.LTIC 
output.   Notice  the  gradual  build  up  of  the  DELTIC  output  as 
time  intervals  of  the  signal  are  introduced.   Note  also  that 
as  time  intervals  continue  to  be  introduced,  the  DELTIC  drops 
out  the  oldest  information  and  represents  only  the  most  re- 
cent time  intervals. 
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As  the  signal  is  "stepped"  into  the  DELTIC,  the  output 
of  the  DELTIC  represents  the  compressed  signal.   Plots  A  to 
F  indicate  the  buildup  of  the  DELTIC.   Plots  G  to  J  show  the 
oldest  information  dropping  out  and  only  the  most  recent  in- 
formation being  represented.   Plot  K  shows  the  complete  time 
history  of  the  signal  and  its  DELTIC  output. 
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VII.  ENERGY  DETECTION 

In  the  discussion  of  random  signals   it  was  noted 
that  in  order  to  apply  statistical  analysis  to  signal 
processing,  the  signals  must  be  "ergodic".   One  property 
of  ergodic  signals  is  statistical  stationarity ,  that  is, 
the  statistical  properties  of  the  signal  are  constant 
in  time.   In  reality,  the  signals  that  must  be  processed 
are  not  strictly  stationary.   The  sea  noise  varies  with 
the  time  of  day,  changes  in  weather,  etc.,  but  the  var- 
iation is  over  a  long  time  period.   This  fact  enables 
us  to  apply  the  statistical  analysis  approach  to  signal 
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processing,  as  xGng  as  ljiC  processing  time  is  Sucr*,  cn^u 
that  the  variation  in  signal  properties  is  negligible. 

Let  us  consider  how  to  apply  statistical  methods  to 
detection  problems.   For  simplicity,  take  the  case  of 
Gaussian  noise,  where  the  probability  distributions  of 

amplitude  for  noise  alone  and  noise  plus  signal  are 

2 
Gaussian  with  equal  variance,  6    ,  but  different  means 

(average  values) . 
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If  we  now  have  a  processor  which  compares  the  amplitude 
of  the  input  to  a  "threshold"  level  set  into  it,  and 
indicates  "signal  present"  if  the  input  is  above  threshold, 
"signal  absent"  if  the  input  is  below  threshold,  then 
we  can  apply  statistical  theory  to  determine  the  oper- 
ating characteristics  of  the  processor. 
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The  area  under  the  S  +  N  curve  above  the  threshold 
is  proportional  to  the  probability  that  the  amplitude  of 
the  input  will  be  above  the  threshold  when  the  signal 
is  present.   Thus,  this  area  is  proportional  to  the 
signal  detection  probability,  P(d). 

P(d) 


The  area  under  the  N  curve  and  above  the  threshold 
is  proportional  to  the  probability  that  the  amplitude  of 
the  input  will  be  above  the  threshold,  even  when  the 
signal  is  absent,  because  of  the  background  noise.   Thus 
this  area  is  proportional  to  the  probability  of  false 
alarm,  P(fa) .  .  .  . 


If  the  threshold  is  varied,  the  probabilities  change 
If  the  threshold  is  raised,  the  probability  of  detection 
decreases,  but  so  does  the  probability  of  false  alarm. 


a 


If  the  threshold  is  lowered,  P(d)  increases,  but 
so  does  P(fa) . 


a 


The  ratio  of  P(d)  to  P(fa)   does  not  remain  constant, 
however.   Plotting  P(d)  versus  P(fa)    as  the  threshold 
varies,  we  obtain  the  following  result. 


P(d) 


P(fa) 


This  is  known  as  a  "Receiver  Operating  Characteristics" 
curve,  or  ROC  curve.   Using  this  curve  we  can  determine 
where  to  set  the  receiver  threshold  to  optimize  P(d) 
and  P(fa)   . 
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As  the  threshold  is  reised,  P(d)  .does  not  change 
much  in  comparison  to  the  rate  at  which  P(fa)    is 
decreasing,  until  the  "knee"  of  the  curve  (between  point  2 
and  point  3)  is  reached. 


P(d) 


P(fa) 

As  we  continue  to  raise  the  threshold  above  this 
value,  P(d)  starts  to  fall  off  rapidly,  while  P(fa) 
changes  very  little.   By  deciding  what  maximum  P(fa) 
can  be  tolerated,  we  can  set  the  threshold  to  yield  the 
maximum  P(d)  possible  for  that  particular  P(fa)    by 
referring  to  the  ROC  curve.   Conversely,  if  a  minimum 
P(d)  is  required,  the  ROC  curve  indicates  the  minimum 
•P(fa)    that  must  be  tolerated. 

Until  now  the  effect  of  changing  the  mean  of  the 
S  +  N  distribution  with  respect  to  the  N  distribution 
would  have  on  P(d)  and  P(fa)    has  not  been  discussed. 
If  the  mean  of  the  S  +  N  distribution  is  raised,  then 
the  ratio  P(d)/P(fa)    is  larger  for  all  threshold  settings 
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The  result  is  shown  on  the  ROC  curve  in  the  following 
diagram. 


P(d) 


P(fa) 


We  define  the  "detection  index"  as  the  quantity 

—      —   22 
(X    -  X  )  /$    ,  which  is  a  measure  of  the  difference 

between  the  distributions.   A  large  d  value  indicates  a 

large  separation  between  the  distributions,  and  thus  a 

curve  with  more  of  a  "knee". 
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P(d) 


P(fa) 


This  diagram  shows  that  if  the  presence  of  a  signal 
drastically  changes  the  input  distribution  mean,  then 
we  can  reduce  the  P(fa)    significantly  while  still 


fj       o    I  o  vrro    T")  C  A 


does  not  drastically  change  the  mean,  large  values  for 
P(d)  are  possible  only  at  the  expense  of  a  large  P(fa)   '. 

We  have  discussed  thresholds,  and  the  way  the  selection 
of  threshold  values  affects  the  probability  of  detection 
and  the  probability  of  false  alarm,  but  have  not  discussed 
how  these  thresholds  are  set.   One  point  that  may  not 
be  readily  apparent  is  that  in  detection  systems  in  which 
the  human  operator  is  involved  in  the  decision  as  to 
presence  or  absence  of  a  signal,  his  perception  of  the 
display  must  be  included  in  computation  of  the  threshold 
value.   Most  systems  in  use  operationally  involve  an 
operator  in  the  decision  process.   The  SQS-23  sonar  is 
an  example.   The  operator  can  vary  the  threshold  to  an 
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extent  by  adjusting  the  brightness  and  intensity  of  the 
display  and  the  gain  control,  but  a  large  consideration 
in  determining  whether  or  not  a  target  is  present  is  his 
ability  to  distinguish  between  target  and  noise  as  pre- 
sented on  the  PPI  scope. 

The  operator  reading  grams  on  an  AQA-7  is  another 
example.   He  can  control  the  gain,  or  the  grey  value  of 
the  background  on  the  paper,  but  his  perception  of  the 
distinction  between  target  and  noise  determines  the 
threshold.   For  this  reason  it  is  difficult  to  determine 
the  ROC  curves  for  systems  with  human  operators,  because 
it  is  difficult  to  determine  the  actual  threshold  setting. 
The  threshold  will  vary  from  day  to  day,  or  even  hour 
to  hour,  with  the  same  operator  and  the  same  hardware 
settings  because  of  the  operator's  state  of  mind,  physical 
condition,  fatigue,  etc.   This  is  what  makes  quantification 
of  the  detection  characteristics  for  a  system  difficult. 
Completely  automatic  systems  are  much  easier  to  evaluate, 
as  one  can  determine  exactly  what  the  settings  are. 

In  the  preceding  discussion  of  energy  detection  and 
thresholding  we  have  seen  how  a  system's  detection  capability 
can  be  determined  from  knowledge  of  the  characteristic 
distributions  of  noise  and  signal  plus  noise.   Implicit 
in  this  treatment  is  the  assumption  that  the  distributions 
of  noise  and  signal  plus  noise  are  exactly  known.   If  these 
distributions  are  exactly  known  we  can  specify  the  per- 
formance of  the  system  for  any  given  threshold  setting 
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in  terms  of  P(d)  and  P(fa)   .   If  we  do  not  know  these 
distributions  exactly,  our  analysis  will  not  exactly 
describe  the  system  performance.   Thus,  the  accuracy  of 
our  prediction  of  system  performance  depends  on  the 
accuracy  with  which  we  can  describe  these  distributions. 
In  designing  real  world  systems  this  is  one  of  the  problems 
facing  the  designer , since  the  distributions  of  noise  and 
signal  plus  noise  are  not  constant  in  time  (over  any 
lengthy  period)  or  with  geographical  position.   For  this 
reason  many  systems  do  not  always  perform  as  well  as 
the  specifications  indicate.   Only  if  the  ambient  condition: 
are  the  same  as  those  for  which  the  system  was  designed 
will  it  perform  as  predicted. 
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VIII.  CORRELATION  DETECTION 

In  the  section  on  random  signals,  power  spectral 
density,  and  noise,  it  was  shown  that  correlation  can 
effectively  improve  the  signal-to-noise  ratio  for  cor- 
related signals  due  to  the  cross  product  terms  in  the 
calculation  of  the  correlation  function.   How  can  this 
information  be  used  to  detect  signals  in  noise?   Let 
us  study  an  active  echo  ranging  system  as  an  example. 

In  this  type  of  system  a  signal  is  transmitted,  and 
then  the  returning  signals  are  processed.   The  range 
of  the  reflecting  object  can  be  determined  from  the  echo 
return  time.   In  a  noise-free  environment ,  there  are 
few  problems  in  this  system,  but  in  most  real  situations 
the  ambient  noise  may  "bury"  some  of  the  returns.   By 
using  correlation,  we  can  improve  the  signal-to-noise 
ratio  and  recover  these  echoes. 

This  system  retains  a  replica  of  the  transmitted 
signal  and  then  inputs  the  received  signal  into  the 
correlator.   The  output  of  the  correlator  is  large  when 
an  echo  is  received,  and  small  when  only  noise  is  present 

What  does  doppler  do  to  the  correlation  of  returning 
echoes?   Let  us  look  at  a  burst  of  CW  signal  transmitted, 
and  the  effect  of  doppler  shifting  on  the  correlation 
function.   If  the  echo  is  not  doppler  shifted,  the  output 
of  the  correlator  will  be  the  autocorrelation  of  the 
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transmitted  signal  at  a  time   ,  which  is  related  to 
the  range  of  the  reflected  object. 
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As  the  returning  echo  is  shifted  in  frequency  by 
doppler,  the  output  of  the  correlator  becomes  the  cross- 
correlation  of  two  closely  related  signals,  but  its  value 
is  never  as  high  as  the  autocorrelation.   The  cross- 
correlation  amplitude  and  the  amount  of  delay  over  which 
it  is  correlated  decrease.  (In  fact,  the  reduction  results 
from  the  decrease  in  the  delay  over  which  the  signal  ib 
correlated  with  its  doppler  shifted  echo.) 
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A  plot  of  the  correlator  output  for  a  number  of  different 
doppler  shifts  on  the  same  diagram,  that  is,  a  T,  Af , 
correlator  amplitude  plot,  shows  the  effects  of  doppler. 
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For  a  given  echo,  this  plot  gives  a  solid  surface 
relating  Af,  t,  and  amplitude. 


If  we  look  down  at  the  T,  Af-plane,  and  plot  the 
contour  line  where  the  amplitude  of  the  function  is 
50  per  cent  of  the  maximum,  we  obtain  the  following 
result . 
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This  diagram  is  commonly  known  as  an  "ambiguity 
diagram",  and  the  reason  for  this  will  become  apparent 
from  examination  of  FM  pulses.   First,  however,  let  us 
consider  a  long  CW  pulse  to  see  how  doppler  affects  the 
correlator  output. 

The  section  on  correlation  and  convolution  shows 
that  the  width  of  the  correlation  function  is  related 
to  the  duration  of  the  time-domain  signal. 
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expected  to  respond  over  a  long  range  of  values  of 
What  about  doppler  and  its  effect  on  a  long  CW  pulse 
being  correlated?   Observation  of  two  signals  at  slightly 
different  frequencies  shows  that  the  longer  the  pulses 
to  be  correlated,  the  less  the  correlation,  since  the 
signals  will  be  almost  correlated  over  only  a  few  cycles, 
but  get  further  and  further  out  of  phase  as  time  progresses 


As  the  difference  in  frequency  between  the  signals 
increases,  the  number  of  cycles  over  which  the  two  are 
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almost  correlated  decreases  rapidly.   Thus  for  long 
CW  pulses,  the  correlator  response  decreases  rapidly 
with  doppler  shifts.   The  ambiguity  diagram  for  a  long 
CW  pulse  is  shown  in  the  following  figure. 


The  length  L  is  proportional  to  the  length  of  the 
CW  pulse.  tR  is  the  delay  corresponding  to  the  actual 
xange  of  the  target. 

For  a  short  CW  pulse,  the  length  of  the  correlation 
peak  50  per  cent  contour  will  be  shorter,  but  the  range 
of  doppler  shifts  over  which  the  signals  remain  essentially 
correlated  is  greater. 
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Thus  the  ambiguity  diagram  for  the  short  CW  pulse 
has  the  following  appearance. 

126 


At 


T+ 


These  two  diagrams  show  why  long  pulse  CW  is  good 
for  searching  and  doppler  determination,  and  short  pulse 
CW  is  good  for  range  resolution.   A  sonar  pulse  of  120 
milliseconds  will  ensonify  a  region  of  water  180  meters 
long.   In  other  words,  the  spatial  length  of  the  pulse 
as  it  travels  through  the  water  ib  130  meters  (1500  meters/ 
second  velocity  x  0.120  seconds  =  180  meters,  the  distance 
the  forward  edge  of  the  pulse  has  traveled  during  the  time 
the  rest  of  the  pulse  was  being  generated) .   The  cor- 
relation peak  50  per  cent  contour,  as  shown  on  the 
ambiguity  diagram,  will  be  proportional  to  the  pulse 
length  T,  and  for  CW  this  turns  out  to  be  a  one-to-one 
ratio.   Thus,  the  correlation  will  be  over  a  period  of 
120  milliseconds,  which  means  that  the  target  range  can 
be  determined  only  to  the  nearest  180  meters. 

In  the  case  of  the  short  pulse  of,  say  5  milliseconds, 
the  pulse  spatial  length  is  7.5  meters,  the  correlation 
peaks  over  5  milliseconds,  and  the  range  resolution  is 
7 . 5  meters . 
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The  effect  of  doppler  on  the  ambiguity  diagram  shows 
that  the  same  type  of  development  can  be  made  to  relate 
frequency  resolution  of  the  two  different  pulse  lengths. 
Because  the  long  pulse  remains  correlated  over  very 
limited  doppler  shifts,  the  frequency  of  the  returning 
echo  can  be  determined  rather  accurately  by  using  a  bank 
of  correlators  comparing  the  echo  to  replicas  with  varying 
amounts  of  doppler  shifts. 
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On  the  other  hand,  use  of  the  same  bank  of  correlators 
with  a  short  pulse  produces  a  response  from  each  of  these 
correlators,  because  the  short  pulse  remains  correlated 
over  a  large  range  of  doppler  shifts. 
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If  we  use  a  pulse  in  which  the  frequency  is  changed 
during  transmission,  which  is  commonly  referred  to  as 
FM  slide,  several  differences  appear. 


FM  "Up-Slide" 
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Because  of  the  frequency  change  during  the  pulse, 
the  returning  echo  correlates  with  the  transmitted  replica 
over  a  short  range  of  time,t,  even  for  long  pulses. 
Thus,  the  ambiguity  diagram  is  very  narrow  in  the  X 
direction.   On  the  other  hand,  if  a  long  pulse  io  used, 
the  Af  range  of  correlation  is  narrow.   One  difference 
between  FM  and  CW  pulses  becomes  apparent  here,  however. 
In  the  case  of  FM  pulses,  a  doppler  shift  will  tend  to 
shift  the  entire  pulse,  in  frequency,  but  the  slide  will 
remain  linear  (provided  it  was  originally  linear) ,  and 
thus  it  will  still  tend  to  correlate  but  at  a  different 
value  of  T.  This  behavior  and  its  effect  on  the  ambiguity 
diagram  are  shown  in  the  next  drawing. 
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Note  that  the  FM  slide  pulse  has  some  of  the  char- 
acteristics of  both  long  and  short  pulse  CW.   At  any 
given  doppler  shift,  the  correlation  peak  itself  is 
narrow,  as  is  the  peak  for  short  CW  pulses.   Thus,  range 
resolution  is  good. 
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Note  also  that  at  a  given  Y  the  range  of  doppler 
shift  over  which  the  signal  is  correlated  is  narrow, 
on  the  order  of  the  range  over  which  the  long  CW  pulse 
is  correlated. 
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The  reason  this  diagram  is  called  an  "ambiguity 
diagram"  can  now:  be  seen.   The  correlator  output,  that  is 
the  f-amplitude  plot,  shows  the  following. 
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Was  the  peak  at  T?  caused  by  a  target  exhibiting 
no  doppler  at  range  equivalent  to  T? ,  or  was  it  caused 
by  a  target  with  5  knots  of  doppler  at  range  T?  -  A  ? 
The  output  looks  the  same  in  either  case.   There  is  a 
"range  ambiguity"  present. 
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One  way  of  resolving  this  is  to  transmit  two  FM  slide 
pulses,  one  with  an  "up-slide",  (increasing  frequency 
with  time),  and  one  with  a  "down- slide" (decreasing  frequency 
with  time). By  then  using  two  correlators,  one  with  an 
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"up-slide"  replica  and  one  with  a  "down- slide"  replica, 
the  ambiguity  can  be  resolved.   The  next  ambiguity  diagram 
shows  the  output  from  both  correlators  for  a  given  range 
target . 
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When  compared  on  a.  T- amplitude  plot,  the  correlator 
outputs  for  a  target  exhibiting  the  doppler  shift  shown 
on  the  ambiguity  diagram  by  the  dotted  line  appear  as 
in  the  following  diagram. 
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By  noting  which  peak  appears  first,  we  can  determine 
if  the  doppler  is  up  or  down.   If  the  "up-slide"  peak 
appears  first,  then  the  returning  echo  has  been  shifted 
down  in  frequency,  and  thus  correlates  at  a  delay,   , 
earlier  than  it  would  if  not  doppler  shifted.   This 
indicates  down  doppler. 
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I£  the  "down-slide"  peak  appears  first,  the  opposite 
is  true,  which  indicates  up  doppler. 
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The  true  range  of  the  target  is  midway  between  the  two 
correlation  peaks,  regardless  of  the  amount  of  doppler 
shift  present.   If  the  target  has  no  doppler,  the  two 
peaks  coincide. 


4* 


A. 

i 


UP«$UD£ 


"T5£*/*A-^LID£ 


The  separation  between  the  up-slide  and  down-slide 
correlation  peaks  is  proportional  to  the  amount  of  doppler 
shift  in  the  returning  echo.   If  a  suitably  calibrated 
display  of  the  two  outputs  is  available,  the  amount  of 
doppler  may  be  read  directly  from  the  presentation. 

Thus,  by  employing  FM  slide  pulses  and  suitable 
processing,  one  can  obtain  the  range  resolution  of  short 
CW  pulses  and  the  doppler  discrimination  of  long  CW  pulses 
from  a  pair  of  long  FM  pulses. 
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Another  type  of  pulse  which  provides  the  "best  of 
both  worlds"  with  respect  to  doppler  and  range  is  known 
as  the  "pseudo- random  noise  pulse",  which  consists  of  a 
"coded"  pulse  with  amplitude  or  frequency  varied  in  a 
pseudo-random  fashion.   Take,  for  example,  the  case  where 
the  amplitude  (or  frequency)  of  the  signal  is  not  a 
simple  sinusoid,  but  a  binary  coded  function. 


P   0 

J.  v. 


^g       w»ir!*S* 


This  pulse  is  correlated  over  a  very  narrow  range 
of  X ,    giving  good  range  resolution,  and  it  is  also  cor- 
related over  a  very  narrow  range  of  doppler  shift  because 
of  its  pseudo-random  character  (Truly  random  signals 
correlate  only  at  a  single  value  of  t.)   The  ambiguity 
diagram  for  a  pseudo- random  pulse  is  shown  in  the  next 
figure. 
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Note  that  because  of  the  correlation  over  narrow 
ranges  of  Af  and  T,  no  range  ambiguity  is  present.   The 
problem  associated  with  this  type  of  pulse,  however, 
is  that  tlic  presence  of  noise  in  the  returning  echo 
causes  marked  degradation  in  the  correlation  of  the 
return,  and  thus  degrades  the  performance  of  real 
systems.   For  this  reason  the  pseudo- random  noise  pulse 
is  useful  in  applications  where  the  signal-to-noise 
ratio  is  high,  such  as  close  range  tracking,  bathymetry, 
etc.,  but  its  usefulness  is  limited  in  long  range  search 
and  tracking,  where  the  signal-to-noise  ratio  is  much 
lower . 

Up  to  this  point  we  have  covered  the  use  of  correlation 
in  active  systems,  and  how  the  use  of  FM  slide  pulses 
combined  with  the  proper  correlation  processor  enables 
us  to  determine  the  amount  of  doppler  in  an  echo  directly. 
Let  us  now  examine  some  applications  of  correlation  in  the 
passive  detection  field. 
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One  difference  between  active  and  passive  correlation 
detection  systems  becomes  apparent  very  quickly.   That 
is,  in  the  active  system  we  compare  the  returning  echo 
with  a  replica  of  the  transmitted  pulse.   In  other  words, 
the  signal  we  are  seeking  is  a  known  quantity  (a  replica 
of  the  transmitted  pulse).   In  a  passive  system,  however, 
no  pulse  is  transmitted,  therefore  we  have  no  replica 
for  comparison  with  the  incoming  signal.   Therefore  we 
cannot  utilize  correlation  in  the  same  manner;  however 
it  is  still  of  use,  particularly  in  passive  direction- 
finding  systems.   The  utility  of  passive  direction-finding 
systems  can  be  seen  by  examination  of  a  typical  system 
and  its  operation.   Consider  a  system  composed  of  two 
hydrophones,  separated  in  space  by  some  distance,  d. 
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If  the  output  of  hydrophone  1  is  used  as  a  replica 
with  which  to  compare  the  output  of  hydrophone  2,  the 
cross-correlation  of  the  two  outputs  can  be  obtained. 
What  does  this  do  for  us?   The  important  property  of  sea 
noise  for  this  system  is  that  it  tends  to  be  isotropic, 
that  is  to  be  the  same  regardless  of  the  direction  from 
which  it  impinges  on  the  hydrophone.   In  addition,  sea 
noise  tends  to  be  uncorrelated  over  distance.   On  the 
other  hand,  most  signals  of  interest  to  us  are  from  a 
single  source,  usually  a  submarine,  and  thus  are  coherent 
to  an  extent,  and  directional  in  nature. 
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Since  sea  noise  is  isotropic,  the  correlator  output 
will  be  small  and  relatively  constant  when  no  signal 
is  present.   Now  consider  what  happens  in  the  presence 
of  a  signal  of  interest.   The  acoustic  energy  traveling 
outward  from  the  source  will  impinge  on  the  hydrophones 
from  some  given  direction. 


Following  a  given  wavefront  as  it  travels  across  the 
hydrophones,  we  see  that  it  will  encounter  one  of  the 
hydrophones  before  the  other  unless  the  direction  of 
impingement  is  perpendicular  to  the  line  joining  the  two 
hydrophones,  or  the  baseline. 
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The  difference  between  the  times  that  the  wavefront 
strikes  the  two  hydrophones  depends  on  the  direction 
from  which  it  arrives  with  respect  to  the  baseline.   If 
the  direction  of  arrival  is  perpendicular  to  this  base- 
line, the  wavefront  strikes  both  hydrophones  at  the  same 
time,  and  thus  there  is  no  time  difference.   If  the 
direction  of  arrival  is  parallel  to  the  baseline,  the 


139 


time  difference  is  a  maximum.   In  fact,  the  time  difference 
is  related  to  the  arrival  direction  by  the  following 
equation: 

At  = 


v  cos  e  ' 

where  v  is  the  speed  of  sound  in  the  medium,  and  6   the 
wavefront  arrival  angle. 


The  cross -correlation  of  the  two  hydrophone  outputs 
for  a  single  wave  passing  the  array,  show?  that  the  delay 
time  7*  at  which  a  correlation  peak  appears  is  a  function 
of  the  wavefront  arrival  time  difference  at  the  hydro- 
phones.  If  the  wavefront  arrival  angle  is  90°,  then 
the  arrival  time  difference  is  zero,  and  thus  the  cor- 
relation peak  will  appear  at  T  =  0.   For  any  other  arrival 
angle,  there  is  a  positive  arrival  time  difference,  At, 
and  the  correlation  peak  appears  at  T  =  At.   This  allows 
us  to  relate  the  correlation  delay  time,  t,  to  the  wave 
arrival  direction  with  respect  to  the  array  baseline. 
There  is  a  problem  of  directional  ambiguity,  as  arrivals 
from  opposite  sides  of  the  array  at  equal  angles  with 
respect  to  the  baseline  produce  the  same  delay. 
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This  ambiguity  can  be  resolved  by  adding  another 
array  to  cross-fix,  or  by  adding  a  third  hydrophone  not 
in  line  with  the  first  two,  to  resolve  the  ambiguity. 


It  should  be  noted  that  although  this  method  works 
well  with  broadband  radiated  signals,  it  has  some  draw- 
backs when  applied  to  signals  of  a  single  frequency  line 
or  those  with  a  very  broad  autocorrelation  function. 
The  problem  is  as  follows.   If  a  signal  consists  of  a 
single  frequency  line,  its  autocorrelation  function  is 
a  sinusoidal  function.   The  question  then  arises  as  to 
which  peak  corresponds  to  the  signal  arrival  direction. 
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For  this  reason,  a  correlation  detector  is  not  well  suited 
for  passive  detection  of  narrow-band  signals,  but  works 
well  for  broadband  signal  detection. 

We  have  now  seen  how  correlation  can  be  of  help  in 
active  systems  and  in  passive  broadband  signal  detection 
systems.   What  can  be  done  about  single- frequency  line 
detection?   In  order  to  answei  Lais  quebLion,  we  examine 
a  system  which  is  similar  to  correlation  in  many  aspects, 
but  is  not  actually  correlation. 
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IX.  BEAM  FORMING 

Beam  forming,  or  electrical  steering  of  arrays,  is 
accomplished  by  a  process  which  is  similar  to  correlation, 
but  with  one  major  dif f erence- -  it  uses  summation,  vice 
multiplication,  of  the  outputs.   We  examine  the  process 
to  see  how  it  operates. 

Let  us  start  with  a  two-hydrophone  array.   As  noted 
in  the  preceding  section,  the  acoustic  wave  arrival  direction 
determines  the  wavefront  arrival  time  difference  at  the 
hydrophones.   With  this  in  mind,  suppose  that  it  is 
desired  to  "steer"  the  array  to  receive  signals  from  a 
given  direction  relative  tc  the  array  baseline  selectively, 
that  is,  to  form  a  "receiving  beam"  at  a  given  orientation. 
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If  we  delay  the  output  from  hydrophone  2  for  a  time 

t,  =  At,  and  add  t,  to  the  time  necessarv  for  the  acoustic 
d  d 

wave  from  the  desired  direction  to  travel  from  hydrophone  2 
to  hydrophone  1,  we  can  selectively  bias  the  array  to 
receive  signals  from  that  direction. 
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Why  does  this  selectively  bias  the  response  of  the 
system  to  signals  from  that  direction?   Consider  a  single 
impulse  arriving  from  the  direction  of  our  selected  beam. 
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Plotting  the  output  of  each  hydrophone  as  a  function 
of  time  produces  the  following  result. 
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By  delaying  the  output  of  hydrophone  2  for  a  time 
td  =  At  and  adding  the  output  of  the  two  hydrophones, 
the  two  impulses  reinforce  each  other  and  the  output 
will  appear  as  in  the  next  drawing. 
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If,  however,  the  impulse  arrived  from  a  different 
direction,  the  time  difference  between  arrival  at  the  two 
hydrophones  will  be  different. 
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If  we  input  the  output  of  hydrophone  2  through  the 
same  delay,  At,  the  plot  of  the  summation  of  the  two 
signals  is  as  shown  in  the  next  illustration. 
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Thus  we  see  that  signals  arriving  from  the  selected 
direction  (time  delay  At)  will  selectively  reinforce 
one  another  in  the  processor  while  those  from  other 
directions  will  not. 

If  we  examine  a  3-hydrophone  linear  array  we  observe 
that  the  effect  is  even  more  pronounced. 
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By  delaying  the  output  of  hydrophone  3  for  a  length 
of  time  2At,  and  the  output  of  hydrophone  2  for  a  length 
of  time  ^t,  a  wave  from  the  desired  direction  will  produce 
the  effect  shown  in  the  next  figure. 
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A  wave  from  a  different  direction  will  give  the 
following  output. 
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It  is  apparent  that  adding  hydrophones  to  the  array 
makes  it  more  responsive  to  signals  from  the  chosen 
direction.   Also,  since  the  output  of  the  system  is  not 
a  correlation  of  the  individual  hydrophone  outputs,  but 
a  simple  summation  of  their  energies,  the  system  will 
work  for  narrow-band  or  single-frequency  line  signals 
as  well  as  for  broadband  signals. 

By  using  a  number  of  hydrophones  in  this  manner  in 
conjunction  with  an  energy  detector,  and  by  setting  the 
threshold  high  enough  to  preclude  response  to  the  un- 
reinforced  signals  from  other  directions,  one  has  a 
very  directional  detection  capability  from  a  series  of 
omnidirectional  hydrophones. 
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